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' Abstract 

<N : 

^ ■ We study the model of a strongly non-linear chain of particles coupled to two heat baths 

^ , at different temperatures. Our main result is the existence and uniqueness of a stationary 

lO I state at all temperatures. This result extends those of Eckmann, Pillet, Rey-Bellet [EPR99a, 

2 ■ EPR99b] to potentials with essentially arbitrary growth at infinity. This extension is possible 

Q>^ . by introducing a stronger version of Hormander's theorem for Kolmogorov equations to vector 

O ' fields with polynomially bounded coefficients on unbounded domains. 

a^ ■ 

Introduction 

-a ; 

Q ■ In this paper, we study the statistical mechanics of a highly non-linear chain of oscillators coupled 
! to two heat reservoirs which are at (arbitrary) different temperatures. We show that such systems 
Q I have, under suitable conditions, a unique stationary state, in which heat flows from the hotter 
\ reservoir to the cooler one. 

These results are an extension of the same statements obtained by Eckmann, Pillet and Rey- 
Bellet in [EPR99a, EPR99b] where it was assumed that the Hamiltonian is essentially "quadratic 
at high energies." Since quadratic Hamiltonians have been discussed much earlier by Lebowitz 
and Spohn [LS77], there is an issue here of whether the quadratic nature of the forces at infinite 
energies is an essential ingredient of existence and uniqueness of the stationary state. Our result 
shows that this is not the case, since we allow for potentials of arbitrary polynomial growth. 
Our models, which are described in Section 1, treat a Hamiltonian of the form 

N 2 ^ 

HsiP, = E (I + + E - Q.-i) , 

i=0 i=l 

describing a chain of particles with nearest-neighbor interaction (Figure 1). This chain is linearly 
coupled to heat baths represented by free fields at temperatures Tj. We proceed then, as in 
[EPR99a], to a reduction to a stochastic differential equation, see (1.2). Associated with it is an 
"effective energy" G, described in (1.6), which is equal to Hg with some quadratic terms from 
the heat baths added. The generator corresponding to the stochastic differential equation above, 
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represented in a space weighted with an exponential of G, will be called K and is the main object 
of study of this paper. 

It is for this generator that we show existence and uniqueness of an invariant state. This will 
be done by first showing that K has compact resolvent (which is really more than needed), and 
then using this result to derive the properties of the invariant measure. Our conditions on Hg are 
spelled out in Section 2 below. They basically say that the coupling between the oscillators must 
be stronger than the single particle potential. This condition might be physically relevant, since it 
implies that transport is favored over storing of energy, but we have not found a counterexample 
when this condition is violated. Furthermore, the interparticle coupling must be convex. 

The main technical insight behind our generalization of the results of [EPR99a, EPR99b] is a 
new, and stronger version of the Hormander theorem for Kolmogorov equations. We will develop 
this in more generality in Section 4, but here we just indicate how we use this result. The operator 
K is of the form 

n 

K^Y.^tX, + X,, (0.1) 
1=1 

where the are smooth vector fields on R'^. For example, see Eq. (2.13), Xq contains terms of 
the foxmp^dq and (5^ y)9p , where V is the interaction. The X^ for z are first order opera- 
tors. Here, dV is polynomially bounded, whereas, in [EPR99a], dV was assumed to be linearly 
bounded. Letting be an adequate inverse power of the effective energy G, one successively 
considers the finite sets of operators 

A_i = {X^, . . . , X„} , = {qoXq, Xi, . . . , Xj^} , 

and then — see Section 5 for the detailed definition — 

Ai = A^^i U [qoXq, Ai^i] . 

We stop this iteration after at most 2N steps, where N is the number of particles in the chain, 
obtaining the set A — .42^+1. We now define the operator A4 as the finite sum 

Al = 1 + X] ■ 

AeA 

This is a generalization to our case of an elliptic operator of the type = 1 — d^, used in 
[Hor85] or A2 = 1 - d^, + J2i used in [EPR99a]. 
With these definitions, one then has the bound 

Proposition 0.1 (Momentum space bound) There is a constant C such that for all f e C^(R'') 
one has in I?: 

l|Ar"/ll < C{\\Kf\\ + 11/11) . (0.2) 

We also derive a similar bound in the conjugate variables: 

Proposition 0.2 (Position space bound) There is a constant C such that for all f e C^{R'^) one 
has in I?: 



WG^fW < + 11/11), 

where e > depends on the asymptotic behavior of the potential V and on N. 



(0.3) 
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Combining these two propositions one easily shows that K has compact resolvent. Then one 
derives from that result the existence of an invariant measure. Its properties are then found adapting 
the techniques of [EPR99a, EPR99b]. 

The remainder of this paper is organized as follows. In Section 1 we describe the physical 
model and in Section 2 we refine the setting and state the results. Section 3 will be devoted to the 
proof of the position space bound (Proposition 2.5). In Section 4, we present in detail the gen- 
eral scheme for studying operators of the form of (0.1), and show the inequality corresponding to 
(0.2). This section is as much as possible self-contained as it presents some independent interest. 
The detailed application of this general scheme to the problem of the chain allows us to prove 
the momentum space bound (Proposition 2.6) in Section 5. In Section 6 we combine these two 
bounds and prove Theorem 2.4 showing that K has compact resolvent and hence discrete spec- 
trum. In Section 7 we show existence, uniqueness, and further properties of the invariant measure 
(Theorem 2.7). 

Appendix A contains some technical estimates used in Section 4. Appendix B contains the 
proof of a result concerning the domains of K and K*. The method used there probably works for 
more general accretive second-order differential operators. Appendix C finally contains the proof 
of a technical result used in Section 7. 
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We will study the model of a (small) classical A^-particle Hamiltonian system coupled to M 
stochastic heat baths proposed in [EPR99a]. The small system without the heat baths is governed 
by a Hamiltonian 



(We stay here with d = 1 dimensional position space for each particle to simplify notation.) The 
heat baths are modeled by classical field theories associated to the wave equation. The fields will 
be called </7j and their conjugate momenta tt^, where the index i ranges from 1 to M. 
The Hamiltonian for one heat bath is given by 



The couplings allowed for the model are linear in the field variables. The total Hamiltonian for our 
model is then given by 



1 The model 



Hs e C~(R'^) . 




M 



H{p, q, Tr,ifi)^Y^ (^HeiTTi, + F^ip, q) I dipi{x)Qi{x) dx^ + Hs{p, q) . (1.1) 
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We assume the initial conditions describe the heat baths at equilibrium at inverse temperatures 
i.e. they are distributed in a sense according to the measure with "weight" 

The paper [EPR99a] explains in detail how, and under which conditions on the coupling func- 
tions Q,i, one can reduce the resulting "big" system to a "small" system, where the heat baths are 
described by a finite number of variables. The price to pay for that is that we are now dealing with 
the following system of stochastic differential equations: 

M 

dqj = d^Hs dt - J2{dpF,)r, dt , j = l,... ,N, 

i=l 

M 

dpj = -d^Hs dt + J2i%^ih ' (1-2) 

i=l 

= -'JiVi dt + XhiFiiP, q) dt - \,^2-fiT,^dw^{t) , i = 1, . . . , M , 

where the are independent Wiener processes. The various constants appearing in (1.2) have the 
following meaning. Tj is the temperature of the heat bath, Aj is the strength of the coupling 
between that heat bath and the small system and 1 /7j is the relaxation time of the heat bath. The 
value of 7j depends on the choice of the coupling function g^^. If we wanted to be more general, we 
would have to introduce for each bath a family of auxiliary variables r,j ^ as is done in [EPR99a]. 
This would only cause notational problems and does not change our argument. 

If we consider a generic n-dimensional system of stochastic differential equations with addi- 
tive noise of the form 

n 

dXi{t) = bi{x{t)) dt + Y^ a^j dwj{t) , (1.3) 
we can associate with it the second-order differential operator JC formally defined by 

-. n n 

>C ^ 2 E + J2 U^) . (1.4) 

i,j=l i=l 

It is a classical result that if the solution of such a system of stochastic differential equations exists, 
the probability density of the solution satisfies the partial differential equation 

dtP{x, t) = (Cp) {x, t) . 

In our case, the differential operator C is given by 

MM M 
i=l i=l i=l 

where the symbol denotes the Hamiltonian vector field associated to the function F. It is 
convenient to introduce the "effective energy" given by 

M 2 

G{p,q,r)^Hs{p,q) + J2{^-Fi{p,q)r,) . (1.6) 

1=1 * 

At this point, we make the following assumption on the asymptotic behavior of G. 



4 The model 

AO. There exist constants dj, C > and ct > 0, as well as constants Cj > 2/Af such that 

Hs{p,q)>C(l + \\pr+\\q\\''), (1.7a) 
F^{p,q)<d^Hs{p,q)+di. (1.7b) 



Remark. This assumption essentially means that the effective energy G grows at infinity at least 
like l + ||r||^ + ||p||"+||g||". This implies the stability of the system, as follows easily from the 
inequality 

|r,F,(p,5)|<.V| + ^5(M, 

s 

which holds for every s > 0. In particular, this implies that exp(— /SG) is integrable for every 
P>0. 

We also define 

M 
i=l 

which is, in some sense that will be clear in a moment, the maximal power the heat baths can pull 
into the chain. We have the following result. 

Proposition 1.1. Assume AO holds. Then the solution ^{t; Xq, w) of (1.2) exists and is continuous 
for allt > with probability 1. Moreover, the mean energy of the system satisfies for all values of 
t and Xq the estimate 

E[G(x(t; Xq, w))] - G(xo) < Wt , (1.8) 

where E[-] denotes the expectation with respect to the M-dimensional Wiener process w. 

Remark. The bound (1.8) allows the energy to grow forever, which would cause the system to 
"explode." But this is not the case for the systems we consider in this paper. Indeed, we will prove 
that the process possesses a unique stationary state. This implies among other features that the 
mean time needed to reach any compact region is finite, and so the energy can not grow forever. 

Proof. A classical result (see e.g. [HasSO, Thm 4.1]) states the following. Assume that the vector 
field 6 of (1.3) is locally Lipshitz and that there exists a confining function G : R" — > R and a 
constant k such that 

{CG){x)<k for all ,t G R" . 

Then there exists a unique stochastic process ^(t) solving (1 .3). The process ^ is regular (i.e. it does 
not blow up in a finite time) and continuous for all t > 0. It satisfies the statistics of a Markovian 
diffusion process with generator C Moreover, we have the estimate 

E[G{x{t;xQ,w))] - G{xq) < kt . 

This result can be applied to our case, if we take for G the effective energy defined in (1.6). An 
explicit computation yields indeed 

M 

jCG^W-J^^iri- miP, q)Y ■ (1-9) 
Moreover, G is confining by AO. This proves the assertion. □ 
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1.1 Definition and simple properties of the semi-group 

In this paper, we will mainly be interested in studying under which assumptions on the chain 
Hamiltonian Hg it is possible to prove the existence of a unique invariant measure for the stochastic 
process ^{t; Xq, w) solving (1.2). Throughout, we will use the notation 

for the extended phase space (p, g, r). This stochastic process defines a semi-group T* on Cq^{X) 
by 

T'f{x,)^E[f{^{t;x„w))]. (1.10) 
This semi-group satisfies the following 

Proposition 1.2. Assume AO holds. Then T* extends to a strongly continuous, quasi-bounded 
semi-group of positivity preserving operators on L^(A'). Its generator L is the closure of the 
operator C with domain Cq°{X). The adjoint L* is the closure of the formal adjoint with 
domain C^{X). 

Proof. The proof will be given in Appendix B. □ 
This in turn defines a dual semi-group (T*)* by 

j {T'f){x)v{dx) = j fix) {{T')*v){dx) . 

The generator of (T*)* is given by the adjoint of C in \3 that will be denoted C^. It is possible to 
check that if the heat baths are all at the same temperature T — we have 

C^Hq — , where A*o(P> ^) — e~^^^^''^''^^ . 

Thus, the generalized Gibbs measure 

d^Q — e~'^'^^''^''^^ dp dq dr — ij,q{p, q, r) dp dq dr , 

is an invariant measure for the Markov process described by (1.2). This confirms our definition 
of G as the effective energy of our system. We want to consider the more interesting case where 
the temperatures of the heat baths are not the same. The idea is to work in a Hilbert space that is 
weighted with a Gibbs measure for some reference temperature. 

We will therefore study an extension Tq of T* acting on an auxiliary weighted Hilbert space 
Hq, given by 

Ho = h^{X, ZQ^e-^f^oG(p,'i,r) dpdqdr) , 

where Zq is a normalization constant and (3^ is a "reference" inverse temperature that we choose 
such that 

l//5o = T^> max{r, | z = 1, . . . , M} . (1.11) 

We have the following 
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Proposition 1.3. Assume AO holds. Then the semi-group T* given by (1.10) extends to a strongly 
continuous quasi-bounded semi-group Tq on Hq. Moreover, TqI — 1 andT^ is positivity preserv- 
ing, i.e. 

Vf>0 if / > . 

Let Lq be the generator o/Tq. Then Lq coincides on C^{X) with C of (1.5) andC^{X) is a core 
for both Lq and Lq. 

Proof. The statement can be proven by simply retracing the proof of Lemma 3.1 in [EPR99a]. 
There are only three points that have to be checked. We define the vector fields h and respectively 
by 

M M 

i=l i=l 
M M 

i=l 1=1 

In order to make the proof of [EPR99a] work, we have to check that 

||div6||(3o < oo , ||div6o||oo < °o ' sup (b + ^bQ)G{x) < oo , 

where b and 6o are considered as first-order differential operators in the last inequality. The diver- 
gence of any Hamiltonian vector field vanishes, and so we have 

M 

I|div6||oo = - J^Ti < oo ■ 

1=1 

The term involving the divergence of 6o can easily be computed to give 

M 

l|div6olloo = /5oZl^i^i 

i=l 

In order to check the last inequality, we compute the expression 

M 

{b + i6o)G(p, g, r) = ^ ^{P,T, - 1) (r, - X^F.ip, q)f . 

i=l 

We see that condition (1.11) on /3o obviously implies /3qT^ — 1 < 0, and so the desired inequality 
holds. 

The domains of Lq and Lq are controlled by the techniques of Appendix B. □ 

We are mainly interested in the case M — 2. The Hamiltonian Hg will describe a chain of A^^-l-l 
strongly anharmonic oscillators coupled to two heat baths at the first and the last particle. In the 
case in which the Hamiltonian Hg can be written as a quadratic function plus some bounded terms, 
the existence and uniqueness of a stationary state for every temperature difference has been proved 
in [EPR99a, EPR99b]. We will extend this result to the case where the potentials grow faster than 
quadratically at infinity. Besides some weak conditions on the derivatives of the one and two-body 
potentials, we will only require that they grow algebraically and that the two-body potentials grow 
asymptotically faster than the one-body potentials, i.e. at large separation the interaction energy 
between neighboring particles grows faster than the one-particle energy. 
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1.2 Notations 

Throughout, the domain of an operator A will be denoted by T>{A). Unless specified, the domain 
of any operator will always be the closure in the graph norm of Co°. For example, if we write 
[A, B], we mean in fact {AB — BA) \ C^, so that the domain of [A, B] can be larger than that of 
Aov B separately. 



2 Setting and results 

In order to set up our model, we need to be able to describe precisely the growth rates of the 
potentials at infinity. This will be achieved with the following function spaces. 

Definition 2.1. Choose a E "R. We call JF^ the set of all C°° functions from to R such that for 
every multi-index k there exists a constant C^. for which 

\\D^f{x)\\<Cj,{l + \\xfY'^ , for all x . 

Definition 2.2. Choose a G R and i G N U {oo}. We call the set of all functions from R" 
to R such that for every multi-index k with \k\ <i, we have f[x) G ^a-\k\- 

Remark. For any a G R, the function 

: R" R 

C2 1) 

X ^ (1 + 



belongs to F^. Moreover, any polynomial of degree n belongs to 



■oo 
n • 



2.1 The chain 




We consider the Hamiltonian 

N 2 ^ 

HsiP, = E ( Y + + E - <l^-l) > (2.2) 

i=0 i=l 
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describing a chain of particles with nearest-neighbor interaction (Figure 1). We sUghtly modify the 
notations used so far. Because there are only two heat baths, we will not use for them the indices 
i E {1,2}, but rather i G {L, R}. Concerning the coupling between the chain and the baths, we 
assume that we can make a dipole approximation, so we set 

Fl^Qo and F^^q^, (2.3) 

in equation (1.1). We will make the assumptions A1-A3 on and V2. 

Al. The potential is in for some n> 1. Moreover, there are constants Cj > such that 

V^{x) > c^P^'^ix) , (2.4a) 
xVl{x) > C2P2"(x) - C3 , (2.4b) 

for all x e R. 

A2. The potential V2 is in some m > n. Moreover, there are constants > such that 

^(x) > c;p2™(x) , (2.5a) 

xV^{x) > c'^P^'^ix) - 4 , (2.5b) 

for all X eR. 
A3. The function 



X 



belongs to .F^ for some £. 

Remark. It is clear that (2.3), together with Al and A2 immediately imply AO. Notice that the 
assumptions Vi G and V2 ^ ^im give bounds not only on the asymptotic behavior of Vi and 
V2, but also of their derivatives. The numbers n, m and £ need not be integers. The generalization 
to a Hamiltonian with Vi, V2 depending also on the number of the particle only creates notational 
problems and is left to the reader. 

An example of potentials that satisfy A1-A3 is 

Vi(x) = - + 2 and V2(x) = (1 + x2)V2 _ cos(x) . 
The effective energy of the system chain+baths is given by 

G{p, q, r)=Hs{p,q) + ^ + ^- - g^yr^ + T , (2.6) 

where we choose the constant F such that G > 1, which is always possible, because n > 1. In 
fact, it is important that the function exp(— /^G) be integrable for any /3 > 0. This could also be 
achieved with for example only one of the one-body potentials non-vanishing, but would cause 
some unimportant notational difficulties. The case n = 1 is marginal, the stability of the system 
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depends on the values of the constants Aj and was treated in [EPR99a]. We will not treat this case, 
but it would not cause any big trouble, as long as G remains confining. 
In the sequel, we will extensively use the notations 

N 

Qi = (li- Qi-i and Q = ^ . 
The system of stochastic differential equations we consider is given by 





Pidt , 


dPo = 


-y/(?o) dt + Vi{q^) dt + rLdt, 


dp J = 


-V{{q^) dt - V^%) dt + V^%^,) dt , 


dpN = 


-Vliq^) dt - V^iq^) dt ^r^dt. 


dVL = 


-iL^L dt + XI'^lQo dt - \l^J2-ilTl dwLit) , 


dvR = 


-IrTr dt + X]aRqN dt - ^27^^71 dw^it) , 



(2.7) 



where i = 1, . . . and j = 1,... ,A^ — 1. Since AO holds, the results of the preceding section 
apply. Therefore, there exists for any initial condition Xo a unique stochastic process ^{t; Xq, w) 
solving (2.7). It obeys the statistics of a Markov diffusion process with generator 

N N /2 g\ 

+ r,d,^ + trO,^ + E (pA. - yi^<l^)^) - E (^j - %j ■ 

i=0 i=l 

We want to prove the existence of a smooth invariant measure with density q, r). It is the 
solution of {T^Yii = 0, where (T*)* is the dual semi-group of T*. To achieve this, we introduce, 
as above, the Hilbert space 

Ho = L2(R2^+4, z-ie-^/^oGb.g.r) dpdqdr) , 

where Zq is a normalization constant and (3^ is a "reference" inverse temperature that we choose 
such that 

l//3o = To > max{r^, Tr} . (2.9) 

Proposition 1.2 holds, so the dynamics of our system is described by a semi-group Tq acting in 
T-Lq with generator Lq, formally given by C. The extended phase space of our system will again be 
denoted hy X = R2^+^. 

For convenience, we would like to work in 7i = L^(A'), so we define the unitary transformation 

(C//)(x) =e^oGW/(x) . 
So Lq is unitarily equivalent to the operator Ly^ : 'D{L'^) H defined by 
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An explicit computation shows that is given by 

where the fonnal expression for the differential operator K is 

N N 



(2.10) 



i=Q 1=1 



Since C^{X) is invariant under the unitary transformation U, it remains a core for both K and K* 
The various constants appearing in (2.10) are given by 



i e {L, R} , 









- (B T 








_ bL bji 




2 2 


a 





ie{L,R} 

We see that condition (2.9) ensures the positivity of the constants a| and a|j, which in turn implies 
that the closure of ReX = (K + K*)/2 is a strictly positive self-adjoint operator. 

The first feature we notice about K is that A3 implies the hypoellipticity of the operators 
K , K* , + K and + K*. We recall that a differential operator L acting on functions in a 
finite-dimensional differentiable manifold A4 is called hypoelliptic if 

sing supp / = sing supp Lf , for all / G V'{M) , 

where V^M) is the space of distributions on Cq^{M). In particular, the eigenfunctions of a 
hypoelliptic operator are C°°. 

The hypoellipticity of the above operators is a consequence of a theorem by Hormander 
[Hor67, Hor85]: given a second-order differential operator 

n 

1=1 

where c : ^ C is a smooth function and the are smooth vector fields. Then a sufficient 
condition for L to be hypoelliptic is that the Lie algebra generated by {Lj | i = 0, . . . , n} has 
maximal rank everywhere. It is not hard to verify that A3 ensures that this condition is verified for 
K, K*, dt + Kanddt + K*. 

Proposition 2.3. If AO and A3 are satisfied, the transition probabilities of the Markov process 
solving (2.7) have a smooth density 

P{t, X, y) e ((0, oo)xX xX) . 
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Proof. This is an immediate consequence of the Kolmogorov equations which state that 

dtP = CP =^ {dt + K- a)U-^P = , 
so C/~^P is an eigenfunction of the operator 9^ + X — a, which is hypoelliptic. □ 

2.2 Main results 

Our main technical result is 

Theorem 2.4. If Assumptions A1-A3 are satisfied, then the operator K defined in (2.10) has 
compact resolvent. 

In order to prepare the proof of Theorem 2.4, we will prove the following two propositions. 

Proposition 2.5. If Assumptions Al and A2 are satisfied, there exist constants C and e > such 
that 

\\G'f\\<Ci\\Kf\\ + \\f\\), for all f eV{K) , (2.11a) 
||G^/||<C(||i^*/ll + 11/11). for all feV{K*). (2.11b) 

Proposition 2.6. If Assumptions A1-A3 are satisfied, there exist constants C, e > 0, a positive 
function Oq : A' — > R and a finite number N of smooth vector fields with bounded coefficients 
such that, for every function f e C^{X), we have 

||A^/||<C^(||i^/|| + 11/11), (2.12) 

where 

N 
i=l 

Moreover, the span the whole o/R^^"*"^ at every point. 

Given Theorem 2.4, we can state and prove the main result of this paper, namely the existence 
and uniqueness of an invariant measure for our Markov process. More precisely, we have the 
following result. 

Theorem 2.7. If Assumptions A1-A3 are satisfied, then the stochastic process ^{t) solving (1.2) 
possesses a unique and strictly positive invariant measure /i. Its density h is C°° and satisfies for 
anypo < t^HPl,Pr}' 

h{x) = /i(x)e-^oGW , 
where h decays at infinity faster than any polynomial. 
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Proof of the bound in position space (Proposition 2.5) 



The above results say that the spectrum of K looks 
roughly like the one schematically depicted in Fig- 
ure 2. We see that it is discrete (compactness of the 
resolvent) and located in the right half of the complex 
plane (m-accretivity). Moreover, it is symmetric along 
the real axis, because X is a differential operator with 
real coefficients. 

Most of the remainder of this paper is devoted to 
the proofs of Theorems 2.4 and 2.7. In the sequel, we 
will always use the notation 



InA 



K 



X, 



Figure 2: Spectrum of K. 



' 



1=1 



where we define 



N N 



i=0 



The operator Xq is antisymmetric, i.e. 



This implies that 



i=l 



Xq — -Xq 



(2.13a) 
(2.13b) 



(2.13c) 



(2.14) 



ItK = J2Xi^i and X^^K-ItK, 



1=1 



(2.15) 



and thus Rei^ is a positive self-adjoint operator. We have one more estimate that will be extensively 
used in the sequel. If / is some function in Cq^{X) and i e {1, . . . , 4} we have 



\\XJf = {f,X:XJ) < {f,^Kf)=Mf,Kf) < WfWWKfW < {\\Kf\\ 
and by a similar argument also 

ii^r/iP<(iwii + "^"^' 



(2.16) 



(2.17) 



3 Proof of the bound in position space (Proposition 2.5) 

First of all, we need a collection of functions belonging to J^q, as defined in Definition 2.1. We 
have the following result. 



Proof of the bound in position space (Proposition 2.5) 
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Proposition 3.1. Let r, p, q and q designate the vectors 

P= {Po,--- ,Pn) > Q= (?!,••• ,Qn) ■ 
Choose a >0 and let /ij. : R*^ ^ R be functions in T^,- Then the functions 

G-'^l''h^{r) , G-'^/''hj,+,{p) , G'-"/(2")/i^^,(g) , and G-'^'^''^'^h^{q) 

belong to T^. 

Proof. We will only sketch the proof of the statement for G~°'/^'^'^^hj^{q). The other expressions 
can easily be treated in a similar way. 

We first notice that G~^{D'^G) is bounded for every multi-index k. This is a straightforward 
consequence of two observations. The first one is that because of the lower bounds (2.4a) and 
(2.5a) of Al and A2 and the expression (2.6) of G, there exists a constant C > for which 

G{p,q,r) > C{r^+p'^ + P^''{q) + P'^"'{q)) , (3.1) 

where P'' was defined in (2.1). The second observation is that, because Vi e J^2n ^2 ^ ^2m-> 
we have for every multi-index k some constant for which 

\D'^G{p,q,r)\<C,{r'+p' + P'^{q) + P'-^{q)) . (3.2) 

Notice that is bounded by a similar argument, in particular because hj^ e JF^. 

We set a — —a/ (2m) and write 

di{G"hj,{q)) = a{G-%G)G<^hj,{q) + G'^d,h^{q) . 

Both terms are bounded by (3.1), (3.2) and the fact that hj^ e T^. It is easy to see that all the 
derivatives can be bounded similarly. The proof of Proposition 3.1 is complete. □ 

Let us define 

The symbol A^ was chosen in order to emphasize the similarity between the proof of Proposi- 
tion 2.5 and the proof of the main result of Section 4, Theorem 4.3. 

Before we start the proof of Proposition 2.5, we notice two more facts. Let us choose a,(3 E R 
with < /5 < 1, and let A, B be two operators of multiplication by positive functions A< B. We 
then have 

{AfAfJ)<{AfBfJ), (3.3) 

as well as the implication 

IIAM/II < C{\\Kf\\ + 11/11) ^ WAfA^fW < C{\\Kf\\ + 11/11) . (3.4) 

Both inequalities are trivial consequences of the fact that A ^ is an operator of multiplication by a 
positive function and the estimate a;^<l + xifa;>0 and s < 1. 
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Proof of the bound in position space (Proposition 2.5) 



3.1 The main tool of the proof 

The main tool in the proof of Proposition 2.5 is the following lemma. 

Lemma 3.2. Let and K be defined as above. Let A and B be multiplication operators repre- 
sented by functions of the form 

h{p,q,r) = ClTl + c^Tr + h{p,q) , h e (R2^+2) . 

Assume moreover that there are exponents and (3^ and positive constants Cj such that the fol- 
lowing estimates are true for every f G C^{X). 

||ArM/|| < CiWKfW + 11/11) , \\^i''Bf\\ < CiWKfW + 

w^r^Afw < C3I1/11 , wAf^Bfw < CM . 

l|Arn^o,A]/ll <^^5(l|i^/ll + 11/11). WAfnXo, B]f\\<Ce{\\Kf\\ + 

If J satisfies the conditions 

1 > «3 + /3i- (3.5) 

/?i + max{/32,/33} „ 
7 > ^ . (3.6) 

7 > min{Q;i + (32, 0^2 + Pi} , (3.7) 

then there exists a constant C such that 

\{[X„B]f,A-'^Af)\<C{\\Kf\\ + \\f\\? > for all feC^{X). (3.8) 

Proof. The proof of this lemma involves some of the commutation techniques developed by 
Hormander [Hor85], but it uses the fact that most operators involved are multiplication operators, 
i.e. they commute. An explicit computation, using (2.6) and (2.13) yields 

[Xo,G]= Tf(^.-^'^i)'' i^i^ G] ^ c^{r J XI -q,), (3.9a) 

ie{K,L} i 

[X2, G] = [X„ G] = 0, [Xs,G]= cj, [r^/Xl - q^) . (3.9b) 

We therefore see that, by Proposition 3.1, we have for i = 0, . . . ,4 

G-^[X,,G]&T^. (3.10) 

Since the X^ are either differentiation operators or multiplicative operators, we have, for any a e 
R, the relation 

G-"[X„G°]=aG-i[X„G']e.Fo, 

and so, since Af = G, 

||Af[X„Ar1|| <oo. (3.11) 

We can now start to bound (3.8). Since [Xg, B] = —X^B — BXq, we can write (3.8) as 

\{[X„B]f,A-^Af)\ < \{BXJ,A-^Af)\ + \{Bf,X,A-''Af)\ 
= T, + T2. 

Both terms will be estimated separately. 
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Term T^. Since we know by (2.15) that Xq — K — RsX, we can write it as 

T, < \{B{^K)f,A^''Af)\ + \{BKf,A^^Af)\ ^ T,,+T,, . 

The term '^^^ be estimated by using (3.7). We indeed have either 7 > + /32, or 7 > + Pi- 
In the former case, we write 

T,, < \\A;'^B\\\\Kf\\\\Ar^'^Af\\ < C{\\Kf\\ + \\f\\r . 

In the latter case, we use the fact that A, B and commute and are self-adjoint to write similarly 

T12 = \{AKf,A-'^Bf)\ < ||Ap^lll|i^/lll|Ar^"^i^/ll < C{\\Kf\\ + \\f\\f . 

Let us now focus on the term T^^. Using the positivity of ReX, it can be written as 

T,, = {{^Ky/^Ar^Bf, i^Ky/^Ar^Af) + ([Ap5, fe/T]/, ApA/) 
= 7i3 + , 

where 

7i,72>0, 7^ + 72 = 7, 

are to be chosen later. We estimate both terms separately. The commutator in T^^^ can be expanded 
to give 

T,, = {Ar'[B,^K]f,Ar'Af) + {[Ar\^K]Bf,Ar'Af) . 

In order to estimate these terms, we recall that ReiT = Yli=i ^t^i- We therefore have 

4 

i=l 

+ {[A-'\Xl]X,Bf,A-'-Af) + (X;[ArSXji?/, ApA/)) 



i=l 



Noticing that [B, X*] is a multiple of the identity operator and that A^ is self-adjoint, we have 

< C{XJ,A-^Af) < \\XJ\\\\A-^Af\\ < C{\\Kf\\ + 



where we used (2.16) and the fact that 7 > 0:2 to get the last inequality. The term T^^^ is bounded 
by C(||ir/|| + 11/11)^ in a similar way. The term t/^^ is written as 

l^f I = \{AnAr\X*]X,f,A-^'ABf) + {Aj^[Ar\X*][X,,B]f,A-^Af)\ 
<C\\XJ\\\\A-^ABf\\+C\\f\\\\A-^Af\\, 

where we used (3.11) and the fact that [X,i, B] is bounded. Now we can bound t/'^^ by Cdl-fC/H + 
11/11)2, using (2.16) to estimate ||XJ|| and (3.7) to estimate ||A^^AB/|| and ||Aj;^A/||. The term 
T^^^^ can be estimated in a similar way. 
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Let us now focus on the term T^g. We can write 



\T,,\<\MKAr^Bf,Ar^Bf)\ 



1/2 



\ i=i 



If we choose 

72 = as , (3.12) 
the terms under the square root are easily estimated by writing them as 

\\X,Ar^Af\\ < ||ApA||||XJ|| + ||[X„Ap]A^^||||ArM/ll 

+ \\^rnXi,A]f\\, 

and estimating the two commutators by (3.11) and (3.9) respectively. 
The term preceding the square root can be written as 

(i^Api?/,Api?/) = (Api?i^/,Ap5/) + ([i^, Api?]/,Api?/) 

= ^15 + ^16 • 

The term T^g can be bounded if we choose 

27i>/3i + /92, (3.13) 

because we have then 

7^15 < \\Kf\\\\Af^B\\\\Af^Bf\\ < C{\\Kf\\ + ||/||)^ . 
In order to estimate the term T^g, we use K = ReK + Xq to write 

Tie = ([Xo, Ap^S]/, ApS/) + ([ReX, Ap^S]/, Ap^S/) 

— 16 "T 16 ■ 

The term T^^ can be estimated by writing it as 

Ti^ = {Arnx„B]f,Ar'Bf) + ([Xo,Ap]A^^ApS/,ApS/) . 
The first term can be bounded by Cdl/T/H + ||/||)^ if we choose 

27i>/3i + /33- (3.14) 

In order to bound the second term, it suffices to have 7^ > Pi, which is the case because of (3.13) 
and the fact that p2> Pi- 

The term T^l^ can be bounded byC(||K/|| + ||/||)^,by treating it in a similar way than the term 
Ti4^. We leave to the reader the verification that no additional conditions on 7^ have to be made. 
This completes the estimate of T^, because (3.12), (3.13) and (3.14) can be satisfied simultaneously 
by (3.6). 
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Term T2. We decompose this term as 

T, < \{Bf,A-^AXJ)\ + \{Bf,A-^[X„A]f)\ + \{Bf, [X„A-'']Af)\ 

= ^21 + ^22 + ^23 • 

Since 7 > ctg + /3i the term 722 is easily estimated by 

T22 < \\A^'^Bf\\\\Arnx„A]f\\ < C{\\Kf\\ + 11/11)2 . 

Noticing that we can assume ai < and /S^ < /Sg, condition (3.7) implies 7 > Since 
[Xq, A^'''] is a function, it commutes with A^^, and so T23 can be estimated writing 

r23<|(A7^^i?/,A?[Xo,Ar^]ArM/)| 

< ||A-'^^i?/||||[Xo,A-^]A7||||ArM/|| < C{\K!\\ ^ 

where we used (3.1 1) to get the last bound. 

We finally bound T^^. Since Xq = K — Re A', it can be expanded as 



The term T^^ can be estimated by writing 

and using (3.7). The term T^x can be written as 

rif = (S/, Ar^(fei^)/) = + (ApS/, [ArM,R2X]/) . 

The term Tyj, has already been estimated. The other term can be treated like the term T^/x- We leave 
to the reader the verification that one can indeed bound it by C ( 1 1 / 1 1 + 1 1 / 1 1 ) ^ without any further 
restriction on 7^ and 73. 

This completes the proof of the lemma. □ 



3.2 The main step of the proof of Proposition 2.5 

By an elementary approximation argument, it is sufficient to prove the inequalities (2.1 1) for / G 
C^[X), since this is a core for both K and K* . Moreover, we will prove only (2.11a). The 
interested reader may verify that the same arguments also apply for (2. 1 lb). 
We want to show that we can find constants e and C such that 

||Af/||<C(||/r/|| + 11/11), for all j ^C^{X) . 

In order to show this, we notice that there is a constant C such that 

N N 

A? < C(l + (r^ - Ai^o)^ + (r^ - XUnY + Y.p" + ^'"(^) + E ^""(^~^)) ^ ^ " 

i=0 i=l 
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The immediate consequence is that 

\\Aifr = {f,Arf)<{f,A'r'Gf). 

It is therefore enough to show that there exists a (small) constant e such that the terms 

iiAr^"(Q)/ii , iiArwii, m-'p-'mh--- 

are bounded by C(||ir/|| + ||/||). 

We are first going to bound the terms involving variables near the boundary of the chain. Then, 
we will proceed by induction towards the middle of the chain. 

The term ||Ar^(ri - \lqo)f\\' We have 

IK^i - Xlqo)fr = l((^L - Aigo)V,/)l < CmK)f, f)\ 

= CmKfJ)\<C\\Kf WW f\\<C{\\Kf\\ + \\f\\)\ 

where we used the fact that a^, 7^ to obtain the first inequality. Since > 1, we thus have the 
estimate 

m-\r^->^ho)f\\<C{\\Kf\\ + \\f\\f 

if we take e <1. 

The term ||Af~Vo/ll- We will prove the estimate 

l|ArVo/ll < + 11/11) > (3.16) 

for Eq <1/ (2m). An explicit computation yields the relation 

[^0, - Ai?o] = h{rL - Xho) - AiPo ■ (3.17) 
Solving (3.17) forpo, we get 

l|Ai°" Po/lP = <AZ'(&L(rL - Aigo) - Xf[Xo,r, - Aigo])/, Af ""W) = - ^ ■ 
The term X^^^ can be estimated as 

l^i'^l < K'lMrL - Aigo)/lll|Ar°-Wll < + WfWr > 

where the last inequality holds because Eq < 1/2. 

In order to estimate Xq^\ we apply Lemma 3.2 with A = Pq and _B = — X^qQ. An explicit 
computation yields [Xq, A] = K/(go) + (o^i) ~ ^l- The term [Xq, B] has already been computed 
in (3.17). Because of Proposition 3.1 and of (3.15), we can choose 

Oil = 1 , Pi = , 

q;3 = 2 - 1/m , /?3 = 1 ■ 

The hypotheses of Lemma 3.2 are thus fulfilled if we choose 7 = 2 — 1/m. We therefore have the 
estimate (3.16) with Eq — 1/ (2m). We have a similar estimate for the symmetric term at the other 
end of the chain. 



Proof of the bound in position space (Proposition 2.5) 



19 



The term II Af ^P'^(qi)f\\. We will prove the estimate 

\\^r"P"'{qi)f\\<C{\\Kf\\ + \\f\\), 

for some £q < Eq. Because of the bound (2.5b) of A2, we can find some constants and C2 such 
that 

<Af°"V-(gO/,/> < Ci|(Af°-V2'(gi)/,gJ>| +C2|<Af°-V,/>| , (3.18) 

2e^ —2 

where we also used (3.3). The second term is easily estimated because A^^ " is bounded if £q < 1. 
We use once again the fact that [Ao,Po] = ^1 (%) + ^'(q^i) ~ to write the first term as 

\{Af°-'V^{q,)f,qj)\ = \{Af°-\[X„p,]-Vl{qo) + r,)f,qj)\ ^ |y « + y/^) + y/^) | . 

The term y/^^ can be written as 

|y(2)| ^ |(Afo-^+^/-l^/(go)/,A-^/™,~J>| < ||Afo-^+^/-l^/(go)/||||A-^/-gJ|| . 

By Proposition 3.1 and the fact that V{ e y^2n-i'> this term is bounded by C||/|p if we take e'q so 
small that 

24<l/n-l/m. (3.19) 
The term y}^^ is bounded similarly by writing 

|y/^)|<||Af°-^+^/-r,/||||Ar^/-g,/||, 

if we impose 

2s'o<l-l/m. (3.20) 

Both conditions can be satisfied because we assumed that 1 < n < m. In order to estimate 

we apply once again Lemma 3.2. This time we have A = qi and B = p^. Using (3.16) and 

Proposition 3.1, we see that we can choose 

= 1/m , Pi^I-Eq, 
a2 = 1/m , = 1 , 

ag = 1 , /Jg = 2 - 1/m . 

By using m > 1, we see that the hypotheses of Lemma 3.2 are fulfilled if (3.19) and (3.20) hold, 
together with e'q < e^jl. Once again, we have the same estimate at the other end of the chain. 

We can now go along the chain by induction. At each step, we go one particle closer towards 
the middle of the chain. We present here only the terms arising when we go from the left to the 
right of the chain. 
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The term ||Af We already treated the case i — 1. Let us therefore assume i > 1. We 

moreover assume that there exist constants > such that we have the estimates 

||Al'-^"Wi/|| < C{\\Kf\\ + 11/11) and ||Af < C(||i^/|| + ||/||) . (3.21) 

We will show that this implies the existence of a constant > such that 

llArWil < ^^(11^/11 + 11/11) • (3.22) 
We use Pi = + [Xq, q^] to write 

The term X^^^ is easily bounded if we write 

1^1'^ I < l|Af^"Wi/||||ArVj|| < C{\\Kf\\ + 11/11)2 , 

where the last inequality is obtained by using Proposition 3.1 and (3.21). We only have to make 

the assumption 26^ < 

In order to estimate the term X-'^\ we apply Lemma 3.2 with A = and B = g,. Explicit 
computation yields [Xq,]?,] = V{{qj) — V2(gj_,_i) — ^2(5^). Using the induction hypothesis (3.21) 
and Proposition 3.1, we see that we can choose 

«i = 1 , = (1 - £-_i)/m , 

Q!3 = 2 - 1/m , /^s = 1 ■ 

If we take < e-_i/(2m), we see that the hypotheses of Lemma 3.2 are satisfied. We thus have 
the desired bound (3.22). 

The term || Af~^P'^(q^^i)/||. We assume that there exist strictly positive constants and £-_i 
such that 

||ArW|| < C{\\Kf\\ + 11/11) and \\A'h-"P"^m\\ < C{\\Kf\\ + ||/||) . 
We will show that this implies the existence of a constant £ ■ > for which 

||Ar'P"^(g,+i)/|| < C{\\Kf\\ + 11/11) . (3.23) 

Expression (3.18) with replaced by q^^^ holds. In order to prove (3.23), it suffices therefore to 
show that 

l(Af^"V2'(gVi)/,?m/)l < C(||i^/|| + 11/11)^ . 

Since, for i > 1 we have [^o^Pj = ^iiQi) ~ ^2 (%+i) ~" ^2(^1)' the preceding term can be written 
as 

\{AT'^-\[X,,p,] + VI{q,) + V^{q,))f,q,^j)\ ^ 1^/^^ + F/^^ + F/'^ | • 
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We impose 2e'^ < 1/n — 1/m. The term y/^^ is then estimated as 

< ||Ar^/"^gm/ll||Af'-'^^^"^W(g.)/|| < c{\\Kf\\ + 

where the last step uses Proposition 3.1 and V{ e ^2n-i- ^ order to estimate the term Y^^\ we 
notice that by the Cauchy-Schwarz inequality and assumption A2, we have 

<C||/||||A;/"^-^P™-ite)Af'->-(gJ/|| 

<c\\f\\\\Ai<-'p"^m\\ 



We can choose 2^^ < so this term can be estimated by the induction hypothesis. The term 
y^^^^ is once again estimated by using Lemma 3.2, this time with A — ^^^^ and B — p^. Using 
Proposition 3.1, it is easy to verify that one can take 



ai 


= 1/m , 


Pi 


= 1-Si, 


^2 


— 1/m , 




= 1 , 


"3 


= 1, 




= 2- 1/m 



It suffices then to choose 2£- < £j to satisfy the assumptions of Lemma 3.2 and get the desired 
estimate. 

It is obvious that this induction also works in the other direction, starting from the other end 
of the chain. It also accommodates to a little bit more complicated topologies, as long as the chain 
does not contain any closed loop. In order to complete the proof of the lemma, we have to estimate 
the last term corresponding to the motion of the center of mass. 

The term || Ai~^P"(Q)/||. Finally, we want to show the estimate 

||Ar^P"(g)/|| < C{\\Kf\\ + 11/11), (3.24) 
for some e. We start with a little computation. We write 

{N + l)qo = Q + (gjv-i - Qn) + 2(giv-2 - Qn-i) + ■■■ + N{qo - Qi) ■ 
Moreover, we have Qi = Qq + {qi — Qo) + ■ ■ ■ + {qi — We can thus write 

This, together with the mean- value theorem, implies the useful relation 

N N 

{N + 1)QVI{Q/{N + 1))=QJ2 Vi{qd + QY, {yi{Q/{N + 1)) - Vi{q,) 

1=0 1=0 

N N N 

1=0 1=0 j=l 



(3.25) 
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where is located somewhere on the Q / {N + 1) and q^. 

In the case of d-dimensional particles, the expression corresponding to (3.25) is 

\{N + 1)QV{{Q/{N + 1))\<{N+1)\Q\\VVM\ 

N N 

+ \Q\J2 sup \V'V,{tQ/{N+l) + {l-t)q,)\J2bij\qj\- 

i=0 *G(0'1) j=l 

The subsequent expressions can be rewritten accordingly. 

We use Al and (3.25) to write the left-hand side of (3.24) as 



II Ar ^p"(g)/ii^ = I {A'r'p'-{Q)f, /) i < c{n + 1) | (Af - v/(Q/(iv + 1)) /, Qf)\ + c\\ff 

N N 

< c 

y(i)+y(2)Vc||/||2. 



iV iV 



The term Y^'^^ can be bounded because Vl' e ^2n-2> so 



N 



k=0 



Thus, Y^"^^ can be split in terms of the form 



The first factor clearly can be bounded by C||/|| if we notice that q i— > -P^"^^(gfc)Q belongs to 
J^2n-i then apply Proposition 3.1. The second factor can also be bounded by C||/|| if we 
impose 

1 1 

< s < — - — , 
2n 2m 

which can be done because we assumed n < m. 

It remains to estimate Y^-^\ We define P = ^^o^'r Since it may easily be verified that 

[^0, P] = Eilo W(%) -r^- r-R, we can write Y^ as 

Y, = (Af -2([Xo, P] + + r«)/, Qf) ^ Y^) + Y^') + Y^') . 

We leave to the reader the verification that the terms Y^'^^ and Y^^^ can be bounded by C||/|p with- 
out introducing any stronger condition on e. The term Y^^^ can be estimated by using Lemma 3.2 
with A = Q and B = P. We have already verified that (3.22) holds for every i, so we can define 

Sp = min{£j I i = 0, . . . , N} . 

This, together with Proposition 3.1, allows us to choose, 

= 1/n , Pi = 1 — Sp , 

02 = 1/n , = 1 , 

= 1 , /^g = 2 - 1/n , 

and thus (3.24) is fulfilled if we choose 2e < Sp. This completes the proof of the lemma. □ 
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4 Generalization of Hormander's theorem 

In a celebrated paper [Hor67], Hormander studied second-order differential operators of the form 

r 

P = Y,L*L^ + L,, (4.1) 

where the are some smooth vector fields acting in R''. He showed that a sufficient condition for 
the operator P to be hypoelliptic is that the Lie algebra generated by {Lq, . . . , L^} has maximal 
rank everywhere. The main step in his proof is to show that there exists a constant £ > and, for 
every compact domain K, C R'^, a constant C)^ such that 

||«||(,) < C^{\\Pu\\ + \\u\\) , V « e Co~(/C) . (4.2) 

In this expression, the norm || • 1 1(g) is the natural norm associated to the Sobolev space H^(R.'^), 



I.e. 



\u\\l,= I \u{k)Wl + kydfk = 11(1 + A)^/2li| 



We base our discussion on the proof presented in [Hor85]. Hormander first defines Qi as the set 
of all properly supported symmetric first-order differential operators q such that for every compact 
domain /C, there exist constants CJ^ and Cjl with 

\\quf < C^F£{Pu,u) + C^'ll^f , u e C^ifC) . (4.3) 

In particular, if we write L* = —Lj + Cj, where Cj is some function, contains all the operators 
of the form 

(L^. - c^./2)A , j>l, 

as well as their linear combinations. It also contains every operator of order 0. Hormander then 
defines Q2 as consisting of the operator {P — P*)/i, as well as all the commutators of the form 
[q, q']/i with g, q' E Qi. For k > 2, he defines Qf, as the set of all commutators [g, q']/i with 
q e Qk-1 and q' e Qk-2- One feature of this construction is that a finite number of steps suffices 
to catch every symmetric first-order differential operator. This is a consequence of the maximal 
rank hypothesis. The main point of Hormander's proof is then the following result. 

Lemma 4.1 (Hormander) Ifq,^ e ande < 2^"^ we have for every /C C R'^ 

Mie-i) < C{\\Pu\\ + \\u\\) , u e Co°°(/C) . (4.4) 

The proof can be found in [Hor85, p. 355]. The result (4.2) then follows almost immediately, 
because the operators id^ all belong to some Q^. Thus there exists some £ > such that 

d 

E < Ck{\\Pu\\ + ll^^ll)^ , u e Co-(/C) , 

which implies (4.2). One of the major problems encountered in this paper is to find a global 
estimate analogous to (4.2), i.e. to find constants C and e such that 

||A^m|| < C{\\Pu\\ + , for all u e C^i'R'^) , 

where A is some modified Laplacean. There are two major difficulties: 
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• If we were to construct the sets Q^. as in [Hor85], they would not necessarily "close" in the 
sense that the successive commutators could blow up, and the whole proof would break down. 
To avoid this we do not necessarily put {P — P*)/i into Q2, but mther Qq^P — P*) /i, where is 
some bounded function. This allows to get decreasing bounds on the successive commutators. 

This problem does not appear in [EPR99a], where the successive commutators are all first- 
order differential operators with constant (or bounded) coefficients. On the other hand, the 
commutator technique is essentially the same as in [EPR99a]. 

• The above construction does not allow to deal with arbitrary symmetric first-order differential 
operators. The reason is that if we want a global equivalent of (4.3), the set Qi is no longer 
allowed to contain products of the Lj and unbounded functions. We thus work with fewer 
operators, which means that we track much more closely the expressions which appear in the 
constructions. 

4.1 General setting 

Let us consider the Hilbert space H = L^(R'^, dx) for some integer d > 1. We define the set €{1-1) 
as the set of closed operators on H and the algebra 03 (H) as the everywhere defined bounded 
operators on H. 

We define T> = C^(R'^), which is dense in H. Let us fix some sub-algebra ^ C ^{H) that 
is closed under conjugation and such that FT> c T> for all F e 5^ (typically ^ is some algebra of 
bounded functions). The advantage of considering C^(R'^) is that every differential operator with 
sufficiently smooth coefficients is closable on it (see [Yos80] for a justification). Moreover, every 
differential operator with smooth coefficients maps T) into itself. This allows us to make a formal 
calculus, i.e. every relationship between operators appearing in this section is supposed to hold on 
T>. The actual operators are then the closures of the operators defined on T>. 

We define £ as the set of all formal expressions of the form 

a^{x)D^ , A; > , a G C°°(R'^) , 

\e\<k 

where denotes the derivative with respect to the multi-index £. By the above remark, any 
element of £ can naturally be identified with a differential operator in €(7Y). 
Consider a differential operator K that can be written as 

n 

x = J]x;x, + Xo, X^e£, j^l,...,n, (4.5) 

1=1 

where Xq is such that 

X*^-X, + g, ged. (4.6) 
We introduce now a definition that will be very useful in the sequel. 

Definition 4.2. Let S G £be a finite set of differential operators and i > a natural number We 
define the set y^{S) as the module on generated by the terms 



S.S^-'-S,, S,eSu{l}, k = l,...,z 
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The elements of 3^5(5) are naturally identified with densely defined closed operators on H. If 
i — 0,we use the convention y^{S) = 

The subscript ^ will be dropped in the sequel when the algebra ^ is clear from the context. We 
construct the sets 

A_,^{X„... ,XJ, Ao = {goX„X„... ,XJ, g,ed, (4.7) 
where the operator Qq is assumed to be self -adjoint, positive and such that 

[9o,Xo\ed. (4.8) 

Let us now construct recursively up to a level R < 00 some finite sets B^, C £,hy the following 
procedure. Assume Ai^i is known. Consider next the set B^^ of all A of the form 

A^ J2 {fsB+Yl fxB[X, 5]) , /b, fxB^d. (4.9) 

We then select a finite subset B^ C Bf'\ The set Ai is then defined as 

Ai = A_,uB,. 

Remark. It is here that our construction differs from similar ones where all elements of Bf^ would 
have been selected. This makes the set of operators which we study much smaller, but then we of 
course have to verify that the operators of interest are really covered by our construction. We will 

make some working hypotheses on the sets Ai. 

HI. The pair {Ar, d) satisfies the following. If A,B e Ar and / e 5^, then 

[A,B]ey\An), A*ey\An), [A/]e^- 



H2. IfAeAi with i > -1, we have A* e y\Ai). 

Remark. Hypothesis HI implies that if X e P(^ij) and F e y'^iAu), then [X,Y] e 
y''~^^~^{Aji). This will be very useful in the sequel. Hypothesis H2 implies that the classes 
y''{Ai) are closed under conjugation. 
We define now the operator by 

A2 = 1+ J2 ■ (4-10) 

This is, in some sense that will immediately be clear from Lemma 4.4, the "biggest" operator con- 
tained in y'^{Aji). The operator A^ is symmetric, densely defined and positive. We will moreover 
assume that 

H3. A^ is essentially self-adjoint on V. 

The powers A" thus exist and are also essentially self-adjoint on P for o; < 2. 
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4.2 Results and a preliminary lemma 

The following theorem is the main result of this section. 

Theorem 4.3. Let K and A be defined as above and assume H1-H3 are satisfied for some R. 
Then there exist some constants C, e > such that for every f eV, we have 

l|AVll<^?(l|i^/ll + 11/11)- (4.11) 
In the sequel, we will write A instead of Aji to simplify the notation. 

In order to prove Theorem 4.3, we need the following lemma, which will be extensively used 
in the sequel. 

Lemma 4.4. Let A, and A be as above and assume HI and H3 hold. If X E 3^|(«4^)' ^hen the 
operators 

Ai^XA^ with /3 + 7 < -j 

are bounded. 

IfY E 2^is such that [F, A^] e y^:g{A), then the operators 

A^[A°^,y]AT with a + P + -f <2- j 

are bounded. 

IfX,Y G £ are such that 

[X, A'] e y^{A) , [r, A2] e y^^{A) and [[A\ X], F] e y^'''-\A) , 

then the operators 

A^[[A",X],r]AT with a + /? + 7<4-j-A; 

are bounded. 

Proof. The proof of this lemma is postponed to Appendix A. □ 

Remark. Lemma 4.4 allows us to count powers in the following sense. Each time we see an 
operator that is a monomial containing fractional powers of A and some operators of y^{A), we 
know that the operator is bounded if its "degree" is less or equal to 0. The rule is that if F G {A) , 
its degree is j and the degree of A" is a. Moreover, every time we encounter a commutator, we 
can lower the degree by one unit. 

Lemma 4.4 also shows that if f e V, A e A and a < 2, expressions such as AA'^f can be 
well defined by 

AA'^f = A'^Af + A"]A-2AV , 

where [A, A"]A~^ is bounded and can therefore be defined on all of H. Similar expressions hold 
to show that any expression of this section can be well defined. 
We are now ready to prove the theorem. 
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4.3 Proof of Theorem 4.3 

The proof uses the commutation techniques developed by Hormander [H6r85] and improved by 
Eckmann, Pillet, Rey-Bellet [EPR99a]. Large parts of this proof are inspired from this latter work. 

Before we start the proof itself, let us make a few computations, the results of which will 
be used repeatedly in the sequel. We first show that we can assume ReX positive. An explicit 
computation, using (4.5) and (4.6), shows that 

n 

= ^X*Xi + ^ , and thus also Xq = K -"S^K + g/2 . (4.12) 

1=1 

Because g E ^, vje can add a sufficiently big constant to Xq to make Re/T positive. This will 
change neither the commutation relations, nor the estimate (4. 1 1). 
Another useful equality is 

goR^K = Re(^o^ + ^i) + ^2 K„ e y\A_i) , (4.13) 

where Ki is a self-adjoint operator such that ^{ggK + K^) is a positive self-adjoint operator. This 
is a consequence of the following two equalities, which are easily verified by inspection 

n 

g,^K = X*g,X^ + K, K, e y\A.^) , 
1=1 

n 

M9oK) = Y,X*g,X^-K, K, e y\A_,) . 
1=1 

We therefore have 

n 

M9oK + K,) = J2x:g,X^. 

i=l 

This proves (4.13). 

Another useful identity will be 

(^0^0)* = -^0^0 + 99o = -90^0 + k> ^0] + 99o 
= -^0^0 + 9o' 9o^^ ' 

where the last equality is a consequence of (4.8). 

We will now verify the estimate (4.1 1) for some vector / e D. In the sequel, the symbol C 
will be used to denote some constant depending only on the operator K. This constant can change 
from one line to the other. We will first prove that A e y^{Ai) with Q <i < R implies 

\\Ay^'^'~^Af\\<Ci\\Kf\\ + \\f\\). (4.15) 

In fact, an immediate consequence of the first part of Lemma 4.4 is that we only have to prove this 
assertion for A e Aj,. The proof will proceed by induction on i. 
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4.3.1 Verification for i = 

We want to verify the estimate 

\\^-'/'Af\\<C{\\Kf\\ + \\f\\), for all AeA- 

The cases A — QqXq and A — Xj with j will be treated separately. 
Tlie case A = Xj. We write 

||A-3/4x./||2 < C\\Xjff < C{f, X*X^f) < C{f, {K + K*- g)f) 

<cMi\Kf)+c\\fr<c\\fmKf\\^ 

This implies the desired estimate. Because X'^ e 3^^(^_i) by hypothesis, this computation imme- 
diately implies the estimates 

||X,/||<C(||X/|| + 11/11), (4.16a) 
ll^;/ll < C{\\Kf\\ + 11/11) , (4.16b) 

which hold for every j > 1. 

Tlie case A — QqXq. We write, using expression (4.12), 

||A-3/M/|p = (^,Xo/,A-3/2A/) 

= {Kf,g,A-^I^Af) + {g,gfA-"'Af)l2 - {{^K)f,g,A-'/'Af) 
= Si + S2 — . 

The terms Si and S2 are easily bounded by C {\\Kf\\ + 1| / 1| )^, using the Cauchy-Schwarz inequality 
and the first part of Lemma 4.4. Using the positivity of Rei^T and the explicit form of K, the term 
can be bounded as 

1^31 = {{^Ky/^f, i^Ky/^goA-'/'Af) 

< \MKf, /) I I {{^K)g,A-y^Af, goA-'/'Af) \ 

, " 1/2 

< V\\Kf\\\\f\\\{ggoA-'/'Af,g,A-^/'Af)/2 + J] WX.g.A-^/^Aff 

i=l 



= VIIA'/llll/ll 



^0 + "^o.* ■ 



i=l 



The term Sq is estimated by simple power counting (the A's contribute for —3 and the A's for 2 in 
the total degree of the expression, hence I^qI < C||/|p). The terms Sq ^ are estimated by writing 

\So,i\ < \\9o^-'^'AXJ\\ + ||[X,,^oA-'/'^]/ll ■ 

The first term is estimated by using (4.16) and power counting. The second term is estimated by 
expanding the commutator as 

[X,,g,A-^/'A] = [X,,g,]A-^/'A + g,[X,,A-^/']A + g,A-^/'[X,,A] , 

and estimating separately the resulting terms. 
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4.3.2 The induction hypothesis 

We shall proceed by induction. Let us fix j > 0, take A e Aj and assume (4.15) holds for i < j. 
Let us moreover define e = 1/4^+^ in order to simplify the notation. Our assumption is therefore 
that 

\\A^-'Bf\\ < C{\\Kf\\ + 11/11) \/Be y\A,_,) . (4.17) 
We will now prove that this assumption implies the desired estimate, i.e. 

||A^- Wll < CiWKfW + 11/11) V A e y\Aj) . (4.18) 

This, together with the preceding paragraph, will imply the estimate (4.15). 

4.3.3 Proof of the main estimate 

Because of the induction hypothesis, we only have to check (4. 18) for A e Aj\Aj_i. By (4.9), we 
can write 

n 

A^ Yl {fBB + fB[9oXo,B] + J2rB[Xi,B]), 

BeAj_i i=l 

with all the / belonging to ^. We have 

n 

W-'Afr^ J2 {{fsB + a9oXo,B] + J2rB[Xi,B])f,A'^-'Af) 

BeAj_j^ i=l 

We are going to bound each term of this sum separately by C(||i^/|| + ||/||)^. 
TermTg. We have 

\Tb\ = |(A2-7sAi-2^A2-iS/,A-U/)| . 

The operators A^^~^/bA^~^^ and A~^A are bounded by Lemma 4.4. Using the induction hypothe- 
sis (4.17), we thus get the bound \Tb\ < C{\\Kf\\ + \\f\\y. 

Term with z ^ 0. We define h = fj^. The term is then written as 

= (B/, X*h*A'--'Af) - {XJ, h*B*A'^-'Af) ^Q.-Q^. 

Term Q^. It can be estimated by writing 

Q, = {Bf, h*A'^-'AX:f) + {Bf, [X*, h*A'^-'A]f) . 

The first term is estimated by rewriting it as 

I {Bf, h*A'^'-^AX*f) I = I {A^'-^Bf, A^-'^'h*A^'-^AX*f) \ 

< \\A^^-'Bf\\\\A'-'^h*A^^-^AX*f\\ < C{\\Kf\\ + WfWf . 
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The last inequality has been obtained by using the induction hypothesis (4. 17), the estimate (4. 16b) 
and the fact that the operator A^~'^^h*A'^^~'^A is bounded by Lemma 4.4. 
The second term is estimated as 

\{Bf, [X:,h*A'^-'A]f)\ = \{A'^-'Bf,A'-'^[X:,h*A'^-'A]f)\ 

< \\A^'-' Bf WW A^-^'[X*,h*A^'-^ A] f\\ . 

The term || A^^^i?/|| is bounded by the induction hypothesis (4.17). The other term can be esti- 
mated by writing the commutator as 

[X*, h*A'^'-^A] = [X*, h*]A'^'-^A + h*[X*, A'^'-^]A + h*A'^'-^[X*, A] . 

The resulting terms are estimated by power counting, using the fact that X* e 3^^ (.4). 
Term Q^. We bound this term as 

\Q,\ = \{XJ,h*A^^'-''AB*f) + {XJ,h*[B*,A'^-'A]f)\ 

< WXJW{Wh*A^'-^AB*fW + Wh*[B*,A^'-^A]fW) 

< ||XJ||(||/i*A2^-2AAi-2=||||A2^-i5*/|| + Wh*[B*,A^'-^A]fW) ■ 

We leave to the reader the not too hard task to verify that it is indeed possible to get the bound 

IQ2I ^ C'dl-^/ll + 11/11)^ by similar estimates as for the term (5]^. 

Term T^. We define h = The term Tg is thus equal to 

n = {[goXo,B]f,h*A^^-'Af) = {g,X,Bf,h*A'^-^Af) - {Bg,XJ,h*A'^-'Af) . 
We use (4.14) to write this as 

TO = - {Bf, h*A'^-'Ag,XJ) + {Bf, g',h* A'^-' Af) 

- {Bf, [g,X„ h*A'^-'A]f) - {Bg,X„ h*A'^-'Af) 
= -U, + U^-U,-U,, 

where gQ G 5^. The term U2 can easily be estimated by 

\U,\ = \{A^^-^Bf,A^-^^g',h*A'^-^Af)\ < \\A^^-'Bf\\\\A^-^^g',h*A'^-^Af\\ 

<C(||i^/|| + 11/11)11/11, 

using the induction hypothesis. In order to estimate the term C/3, we notice that 510^0 ^ 
thus [goXQ, A^] G y'^{A). We can therefore write 

\U,\ = \{A^^-'Bf,A'-'^[g,X„h*A^^-'A]f)\ < \\A^^-'Bf\\\\A'-^^[g,X„h*A'^-'A]f\\ , 

expand the commutator and estimate the resulting terms separately by power counting. We use the 
equality 

Xo^K-ItK + g/2, 

to write the terms Ui and f/4 as 

U, = {Bf, h*A'^-^Ag,{K - (Rei^) + g/2)f) ^ T^^, - T^^, + T^^g , 
f/4 = {Bg,{K - (BcK) + g/2)f, h*A'^-'Af) ^ T„^, - T^^s + T^^^ • 

Each of these terms will now be estimated separately. 
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Terms 3 and T^^q. They are easily bounded like the term U2 by power counting and using the 
induction hypothesis to bound ||A^^~^S/||. In the case of Tg g, we first have to commute B with 
gog/2, but this does not cause any problem. 

Term Tg^i. This term can be estimated by 

The norm of g^A* h?^~'^hJS?~'^^ is bounded by power counting. Using the induction hypothesis 
(4.17), we thus have \Tb^^\ < C{\\Kf\\ + \\f\\f. 

Term T^^. We have the estimate 

IT^^J = \{Kf\glB*h*K^^-^Af) \ < \\K f\\\\g*,B*h*A^^-'Af\\ . 

The second norm can be estimated by writing 

\\g*B*h*A'^'-^Af\\ < \\g*h*A^'-^AA^-^'\\\\A^'-^B*f\\ + \\g*[B*,h*A'^'-^A]f\\ . 

Here, the first term can be bounded by C(|| + ||/||) because, by H2, we have B* e 3^^(w4j_i) 
and so we can use the induction hypothesis. The commutator can be expanded and bounded by 
power counting. 

Term Tb^2' We can write this term as 

T5,2 = {A'^-'hBf, {g,^K)A-'Af) + {A'^-'hBf, [A-' A, g,^K]f) 
= {A^'-'hBf, K,A-^Af) + {A^'-'hBf, [A'^A, go^K]f) 

+ {A'^-'hBf, ^{g,K + K,)A-^Af) 
= Ml + M2 + M3 , 

where the second equality has been obtained using (4.13). These terms can now be estimated 
separately. 

Term M^. We write this term as 

M, = {A^^-^hBf,A-^AKJ) + {A^^-'hBf, [K„A-^A]f) . 
The first term is estimated by using 

K,ey\A_,) ^ l|i^2/ll<C^(l|i^/ll + 11/11) > 

where the implication is a straightforward consequence of (4.16a). The second term can be esti- 
mated by power counting and the induction hypothesis, using the fact that K2 G («4._i), so that 
[K2,A'^A] is bounded. 

Term M2. We use the explicit form of Re/T to write this term as 

M2 = {A'^-'hBf, [A-'A,g,](^K)f) 

n 

+ J2{{^''-'hBf,g,X:[A-'A,X,]f) + {A'^-^hBf,g,[A-^A,XnXJ)) 
+ {A^^-'hBf,g,[A-'A,g]f)/2 

n 

= M20 + + M,2) + . 

i=l 



32 



Generalization of Hormander's theorem 



The term M20 is estimated by using the expUcit form of ReX to decompose it in terms of the form 

\{K^^-^hBf,[A-^A,g,]XtXJ)\ < \\[A-'A,g,]X*m''-'hBf\\\\XJ\\ . 

The norm || [A"M, 5(o]X*|| is finite by Lemma 4.4. The terms \\A'^^^^hBf\\ and \\XJ\\ are bounded 
by C(||A'/|| + 11/11), using the induction hypothesis (4.17) and the estimate (4.16a) respectively. 
The terms M21 and are estimated by power counting and the induction hypothesis. In order to 
estimate the term M^^, we have to commute once more to find 

= {A^^-'hBf,go[A-'A,X,]X*f) + {A^^-'hBf,g,[Xt, [A-'A,X,]]f) . 

The first term is estimated by using (4.16b). The second term is estimated by expanding the double 
commutator and power counting. 

Term Mg. We use the positivity of ^{g^K + K-^) to write 

IM3I = {{M9oK + K,)f'A'^-'hBf, {^{g,K + K,)y/'A-'Af) 

< \Mi9oK + K,)A'^-^hBf,A'^-^hBf)\y'\{Mg,K + K,)A~'Af,A''Af)\y' 

We will now estimate M4, M5 and Mg separately. 

Term M^^. We want to put the operator ggK to the left of /. So we write 

M4 = {A-'hBg,Kf,A'^-'hBf) + {[g,K , A^^-'hB]f , A^^'^hB f) = M,, + M,, . 

The term M^-^ is estimated easily by using the induction hypothesis and the fact that A'^hBg^ is 
bounded. In order to estimate M42, we use the explicit form of K to write 

M42 = {A-^^[goX^,A^'-'hB]f,A^^-^hBf) 

n 

+ Y,{{g,Xt[X,,A^^-'hB]f,A^^-^hBf) + {g^[Xt,A^^-'hB]XJ,A^^-^hBf)) 

i=l 

+ {A-^'[go,A^'-^hB]Kf, A^'-^hBf) 

n 

= M40 + ^(M,3 + M,4) + M^K ■ 

i=l 

The terms M^q and M^^j^ are estimated by expanding the commutator and power counting. The 
term Mj4 can be written as 

\M,,\ = \{A-'^g,[X:,A'^-'hB]X,f,A^^-'hBf)\ 

It is then estimated by power counting, using moreover the induction hypothesis and the estimate 
(4.16). In order to estimate the term M^g, we have to commute once more to write 

M,3 = {A-'^g,[X„A'^''hB]X*f,A^^-'hBf) + {A-^^g,[Xl[X„A^^-%B]]f,A^^-%Bf) . 

The first term is estimated exactly like Mj4. The second term can then be estimated by expanding 
the double commutator and power counting. 
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Term Mg. We write this terni as 

M5 = {A-^hBKJ,A^'-^hBf) + {A-^'[K^, A^'-^hB]f, A^'-^hBf) . 
The first term is estimated using the induction hypothesis and the fact that (4.13) and (4.16) imply 

K,ey\A_,) ^ ||irj||<c(||;^/|| + 11/11). (4.19) 

The other term is estimated by using the fact that [i^^, A^] g y^iA) and [i^^, hB] e y^{A), which 
follows from A_i C A and thus e 3^^^). 

Term Mg. We use the explicit expression for ^{g^K + Ki) to write this term as 

n n 
i=l i=l 

These terms are easily estimated by putting the to the left of /, using (4.16) and estimating the 
commutators. 

Term Tg^g. This is the last term we have to estimate. Using the expression (4. 13) and the positivity 
of Rs(g'o^ + it can be written in the form 

Tb,, = {{M9oK + K,)y/'f, {M9oK + K,)y/'B*h*A^^-^Af) 
+ {KJ, B*h*A^'-^Af) 
= N, + N,. 

These terms are now estimated separately. 

Term N^. We use the Cauchy-Schwarz inequality to write 

\N,\ < \\KJ\\\\B*h*A'^-'Af\\ ^ \\KJ\\\\N,\\ . 

We can estimate A^3 by writing 

B*h*A'^'-^A = h*A'^'-^AA^-^'A^'-^B* + [B*, h*A^'-'^A] , 

and estimating the resulting terms using the induction hypothesis. We already noticed that we have 
the desired estimate for Hi^g/ll • 

Term N-^^. Using the Cauchy-Schwarz inequality, we write it as 

N, < {f,M9oK + K,)fy/^{M9oK + K,)B*h*A'^-^Af,B*h*A'^-^Afy/^ 
< C{\\Kf\\ + \\f\\)\{A-^'{goK + K^)B*h*A^'-^Af,A^'B*h*A^'-^Af)\^/^ 

^ cmn + ii/ii)vi(/i+/2,/3)i < + ii/ii)V(ii/iii + ii/2ii)ii/3ii ■ 

Estimate of 1 1 /g 1 1 . We write it as 

/3 = A2^/i*A2^-2^Ai-4£A4=-i57 + A'^'[B*, h*A'^''^A]f . 

The first term is estimated by using the recurrence hypothesis and the fact that H2 implies B* e 
y^{Aj_i). The second term is estimated by power counting and by using the fact that e < 1/4. 
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Estimate of 1 1 /a 1 1 . We write it as 

/2 = A-^'B*h*A.'^'-^AKJ + A-^'[K^, B*h*K^'-^A]f . 

The first term is estimated using the fact that H-ftTi/H <C'(||i^r/|| + ||/||) and power counting. The 
second term is simply estimated by power counting, and the fact that K^^ e 3^^(^)- 

Estimate of || ||. We use the expUcit form of K to write as 

/i = A-^'B*h*A.^'-^AgoKf + A'^'lg^Xo, B*h*A.^'-^A]f 

n 

+ J2i^-''%X:[X„ B*h*A'^-'A]f + A-'^goX*, B*h*A^^-^A]Xj) 

i=l 

n 

= Qk + Qo + Y.^Qi,i + Qi,2) ■ 

i=l 

These terms will now be estimated separately. 
Term Qj^. We notice that the operator 

A-'^'B*h*A'^'-^Ago 

is bounded by power counting. This yields the desired estimate. 

Term Qq. This term is bounded by C||/|| by power counting, noticing that QqXq e A. 

Term Qi^2- This term can be estimated by power counting if we expand the commutator and use 
the estimate (4.16). 

Term Qi^^. We use once more the trick that consists of putting the X* to the left of /. We write 
therefore 

Qi,i = A-'^5o[^i, B*h*A^^-^A]X:f + [X*[X„ B*h^A^^-^A]] f . 

The first term is estimated by using (4.16b) and expanding the commutator. The second term is 
estimated in a similar way by expanding the double commutator. We don't write the resulting 
terms here, because there are too much of them. They are all bounded by simple power counting 
and by using Lemma 4.4. This completes the proof of estimate (4.18). 

It is now straightforward to prove the theorem. Recall that R is the level up to which the 
are defined. We put e — 1/4^+^, and we write: 

AeA 

= 5^(||A-W||2 + (/,[A%A2-2]^/)) . 
AeA 



The first term in the sum is bounded by using (4.15), the second term by simple power counting. 
This finally completes the proof of Theorem 4.3. □ 



Generalization of Hormander's theorem 



35 



We next note a consequence of this theorem, namely a simple criterion to see if a quadratic 
differential operator has compact resolvent. It is an easy illustration of the technique that will be 
used in the sequel to show that K has compact resolvent. 

4.4 Quadratic differential operators 

Definition 4.5. An operator A : T>{A) — > H is called accretive if it satisfies 

Re(/,^/)>0, for all feV{A). 

An operator A is called quasi accretive if there exists A e R such that A + X is accretive. It is 
called strictly accretive if there exists A > such that A — Xis still accretive. 

If — y4 is accretive, A is called dissipative. An operator A is called m-accretive if it is accre- 
tive and if {A + A)~^ exists for all A > and satisfies \\{A + A)~^|| < A~^. The expressions 
m-dissipative, quasi dissipative, etc. are defined similarly in an obvious way. An equivalent char- 
acterization of m-accretive operators is that they are accretive with no proper accretive extension. 

It is a classical result (see e.g. [Dav80]) that the quasi m-dissipative operators are precisely the 
generators of quasi-bounded semi-groups. An immediate consequence is that if an operator A is 
(quasi) m-accretive (m-dissipative), its adjoint A* is also (quasi) m-accretive (m-dissipative). 

Proposition 4.6. Let Ti. be a Hilbert space and C be a dense subset ofTC. Let K : T>{K) ^ H be 
a quasi m-accretive (or quasi m-dissipative) operator and let : V^A?) Tibe a self-adjoint 
positive operator such that C C V{A?). Assume moreover that C is a core for K, that A? has 
compact resolvent and that there are constants C > and < £ < 2 such that 

||A7||<C(||X/|| + 11/11), for all fee. (4.20) 

Then K has compact resolvent too. 

Proof. By assumption, there exists a constant A > such that -I- A is strictly m-accretive. 
Moreover, (4.20) with K replaced hy K + \ holds if we change the constant C. Since C is a core 
for K, a simple approximation argument shows that ^{K) C 'D{K^) and that (4.20) holds for 
every f eV{K). 

This immediately implies that {K -\- \)*{K -\- A) has compact resolvent. Since {K -\- A) is 
strictly m-accretive, it is invertible and the operator 

{{K + XY{K + A))-^ = {K + \)-\{K + X)-y , 

is compact. Moreover, we know that {K-\-\)~^\s closed, so we can make the polar decomposition 

{K + X)-^ = PJ , 



with P self -adjoint and J unitary. Thus is compact. By the spectral theorem and the charac- 
terization of compact operators, this immediately implies P compact, and thus also PJ compact. 
Thus K has compact resolvent. □ 
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We now consider V, — L^(R'^) and — {\I \ A e R}, where / is the identity operator in H. 
We define the formal expressions 

— (Xi, . . . , x^) , 

Let ^4 : R'' ^ R'' be a linear map and 

B = {6^ e R'' I i = 1, . . . , s} , C = e R^ I i = 1, . . . , , 

two vector families. Let us consider the differential operator K defined as the closure on (R'^) 
of 

s t 

K = -Y^ d^b.bjd^ + x'^cjcjx + x'^Ad^ . (4.21) 
i=i j=i 

We are interested in giving a geometrical condition on A, B and C that implies the compactness of 
the resolvent of K, and therefore the discreteness of its spectrum. It is possible to prove that K is 
quasi m-accretive. Just follow the proof of Proposition B.3, replacing G{x) by x'^x. 
We have the following result. 

Proposition 4.7. A sufficient condition for the resolvent of the operator K defined in (4.21) to be 
compact is that the vector families 

[j (A^)^B and [j A^C (4.22) 

Ar>o iv>o 

span the whole space R". 

Remark. The intuitive meaning of this theorem is that we can apply Hormander's criterion in both 
direct and Fourier space to obtain an estimate of the form 

Ili/VII < + 11/11) , H^-d^d, + x^x. (4.23) 

It is well known that H has compact resolvent. By Proposition 4.6, (4.23) implies that K has 
compact resolvent. 

Proof. We have the following relations 

[x'^Ad^,b'^d^] = Y^[xiaijd^,,bf,d^J = ^ ftj^^, ^^Ja^j-^^. 

i,j,k i,j,k 

= - h^kiaijdx. = -b'^Ad^ , 
[x'^Ad^,c'^x] = ^[a^jOj^-Sj.^, c^x^] = y^^Xidijld^^, Xf.]cf^ 

— 5^ ^i'^ij^jk^k — C^-^^X . 
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We take — 1, so we have Aq — A_i U {x^Ad^}. We construct the remaining Ai by 

B, = [x^Ad,,B,_,]. 

It is very easy to verify HI and H2, because the assumptions we made on A, B and C imply that 
y^{A) contains every operator of the form b^d^ or c^x. We have moreover 

(b'^dj* - -b^d^ and (c^x)* = c^x . 

It is well-known that H3 concerning the essential self-adjointness of the constructed in Theo- 
rem 4.3 holds. Finally, it is straightforward that satisfies >CH, where H is the "harmonic 
oscillator" defined in (4.23). 

This proves the validity of (4.23), and hence of the assertion. □ 

The interested reader may verify that Proposition 4.7 is quite stable under perturbations. A 
similar result indeed still holds when the coefficients b^ and Cj are not constants, but functions in 
Tq. This is precisely what was proved in [EPR99a]. 

5 Proof of the bound in momentum space (Proposition 2.6) 

This proposition is an application of Theorem 4.3. It is just a little bit cumbersome to verify the 
hypotheses of the theorem. In this section, the symbol K will again denote the operator defined in 
(2.10). 

We choose 

5^ = -^0 ' 

which is simply the set of bounded smooth functions with all their derivatives bounded. It is trivial 
to check that is an algebra of closed operators. Moreover, they are all self-adjoint. We also define 

V = C^{X). 

In this section, we will first construct a set A according to the rules explained in Section 4. 
Then we will check that H1-H3 are indeed satisfied, so we will be able to apply Theorem 4.3. 
This will prove Proposition 2.6 almost immediately. 

Before we start this program, we write down once again the definition of Xq, as it will be used 
repeatedly throughout this section: 

N N 
i=0 i=l 

5.1 Definition of A. 

We choose an exponent a < —3/2 — £/ (2m) and we let Qq be the operator of multiplication by 
G". It is clear that Qq is self-adjoint and positive. Moreover, we recall that 

[Xo,G"] = aG"G-i[Xo,G]e.Fo, 
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and so we have [Xq, g^] e ^. The set Aq is defined as 

A = {cl\, CrO,.^, G^Xq} U a , with A = {ai(ri - Xlq^), a^^TR - XI^n)] ■ 
Before we define the sets A, we need a few functions. Let i > be a natural number. The 



functions V^^ and V^' are defined respectively by 



Vr{q)^Vi'mV^{q,-,)-...-Vi'{q,), 
V^\q) = l^2"(?iv+i-0 • ■ ■ ■ • ^2"(giv-i)^2"(?iv) 



It is useful to notice that 



d,yi:\q) 



0, 

{Vi"{q,)Vi'{q,)-')V^\q), 

I [VrmVi'm-' - Vi'\q,^,)Vi\q,^,)-')V^\q) , 

.{{) 



if j > i , 

if j = , 

if j = i , 

otherwise . 



(5.1) 



There are symmetric relations for the derivatives of \ At this point, we use assumption A3 to 
write 



(5.2) 



This implies 



N 



^ G''Y.PJf^J^L^^ = f\vi!\q) , f^e:s, (5.3) 

j=0 

because of Proposition 3.1 and by the choice a < —3/2 — £/ (2m). Moreover, we notice that 

Still because of Proposition 3.1. One more thing we have to remember is (3.10), which implies for 
example that there exists a function f^E^ such that 

[CXq, G^] = , for any /? e R . 

We are now ready to complete the construction of A. 

5.1.1 Definition of A.^ and A2 

We verify that in the case of our model, we can find functions and fxn (4-9) such that 

AMo - {G^'d^,, G-d^J , 
AMi = {G"^^,, G^'^d^J . 

Considering the elements of Ai, we see that it is indeed possible to write 

G^d^^ = cl'[G-X„ c^d^J - G-\d,^G)G-X, + G^h^d,^ , 
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and a similar relation concerning C^dp^. The operators G~^{d^^G) and C^b^cl^ belong to ^, so 
we succeeded to construct Ai according to (4.9). 

Let us now focus on the elements of A2. We can write 

G^'^d^^ = [G"Xo, G'^dJ - G'^-\dpG)G'^X, - af.G^dp^ , 

and an equivalent expression at the other end of the chain. Since G°-~^{dpjG) e ^ and /q G 5^, we 
succeeded to construct A2 according to (4.9). 

5.1.2 Definition of A^i_-i and 
For i > 1, these sets are defined by 

A.-i\A.-2 = {G^'-^'^'^V^'^dp^, G^'^^'^'^V^^dp^J , 

We repeat this construction until i = A?^ — 1, i.e. we do not stop at the middle of the chain, but we 
go on until we reach the other end. We want to check that these sets were constructed according to 
(4.9). In fact, we will see that any element A of Aj\Aj^i with j >2 can be written as 

A^[G''X^,B]+D , BeAj_,, Dey\Aj_,). (5.4) 

We will verify this only for 2 < i < — 2. We let the reader verify that (5.4) is also valid for the 
remaining sets. 

Let us first take j = 2i - 1 and ^ = G^^'+^^^'^V^"'^ dp,. We choose B = G'^^'^V^"'^^ d^,_^ e Aj_i 
and write 

[G"Xo,i?] = h-yG^^'^vt'^ d^^_^ - G'^^vt'^G-\d^^_G)G-X, 
+ 2tafoB + G('^+''>-vt'\x,,d^J. 

The first three terms belong to y^{A2i-2) and can thus be absorbed into D. The last term can be 
written as 

G('^+'^<^vt'\x„ d^J = G^"(V^i"(%-i) + V^'iq,) + V^{q,_,))G('^-'^-vt'^dp^_^ 

+ G'-V^'iq)V^{q)G('^-'^-vt'^dp^_^ 
+ GC'^+'^'^V^^dp^ . 

The first two terms also belong to iy^(^2i-2)> so they can be absorbed into D as well. The remain- 
ing term is 

Q{2i+l)ay^)Q^^ = A , 

thus we have verified that A can be written as in (5.4). The procedure to get the symmetric term 
from the other end of the chain is similar. 

We take now j = 2i and A = G^'^'+'^^'^V^^ d^,. We choose B = G^'^'+^^'^V^^ dp, e J^H-^j-i) 
and write 

+ {2i + l)a/o5 + G'(2^+2)«y^*)a^^ . 

The first three terms belong to 3^^(^2i 1) ^i^d can be absorbed into D, so we verified that every 
element of A can indeed be written as in (4.9). 
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5.2 Verification of the hypotheses and proof 

In order to be able to apply Theorem 4.3, we verify the hypotheses H1-H3. 

Verification of H2. We want to check that A e Aj implies A* e y^{Aj). By Proposition 3.1, we 
can easily verify that A\A C Cq. But we know that 

A e Cq^ A* ^ -A + g , 

and so H2 holds for ^\^. The elements of A being self-adjoint, H2 holds trivially. 
Verification of H3. The operator can be written as 

A'^-J2d,a,,{x)d, + V{x). 

It is well-known that if a^j and V are sufficiently nice, such operators are essentially self-adjoint 

on C^{X) (see e.g. [Agm82, Thm. 3.2]). 

Verification of HI. Let us define S1,q C £, as the set of first-order differential operators with 
coefficients in JFq. We first verify that 

AeA,fed ^ [A/]e^- 

This is trivial, noticing that .A C £o U ^ and [£o, d] = [A,d] = 
We now verify that 

AeA =^ A* e y\A) . 

This is also trivial, because A e 2q => A* — —A + g, with g e J^q. Moreover, the elements of A 
are self-adjoint. 

Finally, we want to verify that 

A,BeA =^ [A,B]ey\A). 

This is a little bit longer to verify. 

Concerning the commutators of the elements of A with the other elements of A, the statement 
follows easily from the fact that if F : R" ^ it! is linear and A e 2q, then [A, F] e J^o = ^. 
Moreover, the commutator between two multiplication operators vanishes. 

Concerning the commutators between the dj. and the other elements, we notice that they com- 
mute with the functions V^'^ (g) and V^^ {q). Moreover, we have for example 

[a,^,G^]=7(G-M\,G])G"^, if 7eR, 

and [dj.^ , G] belongs to ^. It is straightforward to verify that this implies the desired statement. 

Concerning the commutators of G^'Xq with the other elements of A, the statement has already 
been verified by the construction of A for every operator, but those in A2n-2\^2n-3- These 
operators are of the form 

A = , 
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and a similar term at the other end of the chain. We can make a computation very similar to the 
one we made when we constructed ^2i-i' to show that 

[G"Xo, A] = G^^^+^^^l^i'^^o'p^ + C, Cey\A) . 

But e 5, so [CXo, A] e y\A). it remains therefore only to verify the statement for 

commutators between elements of A\Ao. We can divide these commutators in three classes. 

Both operators contain a dp. We notice that these operators can all be written in the form 
G'"<Wj(g)9p . The commutator between two such elements is given by 

- G-\dj,G)G''^Wj{q)G^iW,{q)dp, . 



J 

I 



Both terms belong to y^{A), because G-^{dpG) e d- 

One operator contains a dp, one contains a d^. Let us compute the commutator between 

Q{2i+2)aY{i)Q G'(2^+i)°T4''^a„ . We have 

^Q{2i+2)ay(i)Q^^^ G(^j+^)°'V^^'^ d^^] = G('''+''^°'Vi!\dgG)G-^G'^''^+^Wif^dp 

+ G(^'+^^"Vj^\dpG)G-^G^^'+^^''Vl!^d, 
+ G^'^'V^^G^'^fi-G^^^+^^^V^'^dp, . 

All those terms belong to y^{A). The computation is similar if we take for example 
instead of . 

Both operators contain a d^. The computation is similar to the preceding case and is left to the 
reader. 

It is now easy to give the 

Proof of Proposition 2.6. We have just verified that the hypotheses of Theorem 4.3 are satisfied. 
We apply it, so we have the estimate 

IIAVII < C{\\Kf\\ + 

where A is given by 

A = 1 + ^ >1M . 

AeA 



It is easy to see that A has exactly the form (2. 12). This completes the proof of Proposition 2.6. □ 
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6 Proof of Theorem 2.4 

It is now possible to prove that the operator K has compact resolvent, which is one of the main 
results of this paper. Before we start the proof itself, we need two preliminary results. The first 
one states 

Lemma 6.1. Let A be the closure in L^(R") of the operator acting on C^(R") as 

N 
i=l 

where the are smooth vector fields with bounded coefficients spanning R" at every point and 
is a smooth positive function. 

Let V : ^ be a continuous function such that for every constant C > 0, there exists 
a compact Kq C R" with the property that V{x) > C for every x e W\Kfj. We moreover 
assume that V{x) > 1. Define the operator H as the closure in L^(R") of the operator acting on 
f e Co°°(R") as 

{Hf){x) = {Af){x) + V(x)f(x). 

Then the operator H is self-adjoint. 

Suppose V and the are such that the function 

N 

2aoy + 5]((L* + Lj[L„y] - [L,,[L,,V]]) (6.1) 

i=l 

is bounded. We then have the estimate 

if, H^f) < if, A^f) + (/, V^f) + C{f, H^-'f) , < £ < 1 , (6.2) 
which holds for any f e C^(R"). 

Proof. The result concerning the self-adjointness of H and of A is classical, we will not prove it 
here. The interested reader can find a proof in [Agm82, Thm. 3.2]. 

We use the fact that if T is a strictly positive self-adjoint operator and q; = 1 — £e(0,l),we 
can write 

T-" = C„r.-«(. + T)-d., C^='^. 

Jo TT 



and thus 



T'^C„, I z'-' -dz. 

z + T 



a 



Moreover, a core of T is again a core of T*^, so (6.2) makes sense. For a proof of these statements, 
see [KatSO, §V.3]. This allows us to write inequality (6.2) as 



oo 







\ z + H / Jo \ z + A / Jo \ z + V / 

(6.3) 



foo I 1 

v£-l/ " 
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In order to prove (6.3), let us first show that the operator AV + VA is lower bounded. This is an 
immediate consequence of (6.1) and the equality 

L*L,V + VL*L, = 2L*VL, + (L, + L*)[L„ V] - [L„ [L„ V]] , 

which is easily verified, using the fact that + L* is simply a function. Therefore, there exists a 
constant C > such that 

{g, {AV + VA)g) + C{g, g) > , e C^{R-) . 

Since if > 1 in the sense of quadratic forms, we find 

{g, (AV + VA)g) + C{g, (z + H)g) > , 

which holds for every 2; > 0. Since A and V are positive self -adjoint operators, this immediately 
implies 

(5, ^J^^a) + (5, ^^^^9) + {9: i^V + VA)g) + C{g, {z + H)g) > . (6.4) 

We can easily check the following identities 

VA{z + A)-^V + AV = {z + H)A{z + A)-^V , 
AV{z + V)-^A + VA = {z + H)A{z + V)-^A . 

Inserting this in (6.4), we get 

{g, {z + H) {A{z + A)- V + V{z + V)-^A)g) + C{g, {z + H)g) > , 

and thus 

{g, {z + H)Hg) < {g, {z + H){H + A{z + A)- V + V{z + V)-^A)g) + C{g, (z + H)g) . 

We can check the equalities 

A{z + A)-^V + A = A{z + A)-\z + H) , 
V{z + V)-^A + V ^V{z + V)-^{z + H) , 

which allow us to write 

{g, {z + H)Hg) < {{z + H)g, {A{z + A)-i + V{z + V)-^){z + H)g) + C{g, {z + E)g) . 
Let us define / = {z + H)g. This immediately yields the estimate 




(6.5) 



which holds for any / in W = {z+H)C^ (R") . But we know that (R") is a core for H, therefore 
W is dense in L^(R"). Since the operators appearing in (6.5) are all bounded, the inequality (6.5) 
holds for every / e L?(R") and thus in particular also for / e Co°(R"^). This implies the wanted 
estimate (6.3). □ 
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The second result we want to use is 
Proposition 6.2. Let A, V and H be as in Lemma 6.L Then H has compact resolvent. 
Proof. We know that A is a positive self-adjoint operator, so 

T= (A + l)-i 

exists and ||T|| < 1. The proof of compactness is a modification of the standard proof of the same 
theorem with A replaced by the true Laplacean A, which can be found e.g. in [Agm82]. It is 
based on the fact that if x is a function with compact support, then the multiplication operator x is 
relatively compact with respect to A. We want to prove that xT is a compact operator, i.e. that the 
closure of 

Y^{xTf\feC^{^^) and ||/|| < 1} 

is compact. 

Let us define K, = suppx- By hypothesis, K, is compact. Moreover, we have Y c Cq^{IC). It 
is well-known that if /C is a compact domain of R", then the set 

{u e C^(/C) I \\u\\ < 1; {u,Au) < 1} 

is compact (see e.g. [RS80, Thm. Xin.73]). This implies that Y is compact if we are able to prove 
that there are strictly positive constants e, and C2 such that u eY implies 

< Ci and {u,Au) < C2 . 

We take any element u in y and write it as w = X^f- We have 

M< \\x\\oo\\T\\ 11/11 <ci. 

Recall that we assumed the vector fields L,j appearing in the construction of A span R" at any point 
and that Oq is a strictly positive function. Together with the compactness of the support of u, this 
implies that there are constants C and such that 

\{u,Au)\<C\{u,Au)\=C\{u,AxTf)\ < C\\u\\\\AxTf\\ 

< CllxAT/ll + C||[A,x]r/|| < k, + C\\[A,x]Tf\\ , 

where the last inequality is a consequence of T = (1 + A)~^. We therefore only need to bound the 
term containing the commutator of A and x- Explicit calculation yields 

N N 

X] = E x]Li + [L„ [L„ xW + m + L,) [L„ x]) = J] + % , 

i=i 1=1 

where the rf^ are bounded functions with supp rf^ C /C. So the only terms that remain to be bounded 
are of the form ||77jLjT/||. As 77^ is bounded, it is enough to bound ||LjT/||. We have 

||L,T/||2 = (T/, L*L,Tf) < {Tf, AT/) < \\f\\' . (6.7) 
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This completes the proof of the statement about the relative compactness of x- 

This implies that we can add to H any function with compact support without changing its 
essential spectrum (see [RS80, Thm. XIII. 14]). But the assumption we made concerning V and 
the positivity of A imply that for any constant C, we can raise the spectrum of if + % above C by 
taking for % a smooth function satisfying 

C X E ICq , 
d{x,]Cc)>l. 

Therefore, the essential spectrum of H is empty and thus H has compact resolvent. □ 



It is now easy to give the 

Proof of Theorem 2.4. By Proposition 2.5 and 2.6, we can choose a constant e small enough to 

have, for every / G C^{X), the estimate 

l|A^/ll<C(l|i^/ll + 11/11) and ||G7||<C(||ii'/|| + 11/11). 
We moreover define 

H = A + G. (6.8) 

By the proof of Proposition 2.6, we see that the assumptions of Lemma 6.1 are satisfied. We can 
thus write 

WH^W = (/, H^'f) = (/, (A + Gf^f) < (/, A^^) + (/, G'^f) + C\\fr 
< II A^/IP + IIG^/IP + CII/IP < C(||X/|| + 11/11)2 . 

Because G is confining, we can apply Proposition 6.2 to see that H, and therefore also H^, have 
compact resolvent. Therefore Corollary 4.6 applies, showing that K has compact resolvent. □ 

Remark 1. The proof still works under slightly weaker assumptions. The coupling between the 
ends of the chain and the heat baths does not have to be of the dipolar type. It is enough for 
example that and belong to some .F^ with (3 < n. Moreover, the potentials and V2 can 
be different for each particle. We only have to impose that assumptions A1-A3 can be satisfied for 
every particle with the same constants £, m and n. 

Remark 2. Throughout this paper, we restricted ourselves to the one-dimensional case, i.e. each 
particle had only one degree of freedom. It is not very hard to generalize the results of this paper 
to the d-dimensional case. It is straightforward to generalize assumptions Al and A2, where V is 
now a vector. In assumption A3, the inverse of V2 has to be read as the inverse matrix. A matrix 
or vector- valued function is said to belong to if each of its components belong to Tp. 

The only point that could cause some trouble is the expression (5.1), because the ^"(gj) are 

now matrices which do not commute, so the expression for d^ V^^^ (q) will show terms of the form 

■■■■■ Vi"%) Vi'{q,) , 
where V2' is a trilinear form. Such a term can we written as 

Vi'mVi\qr-i) ■■■■■ Vr{q,)V^%^,)-' ■ . . . ■ V^iA^'vl^ . 

If we want to get expressions similar to (5.2) and (5.3), we have to make jaj very big (of the order 
of A^"), but this is not a problem. 
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Remark 3. One important assumption was that m > n,m other words, the interparticle coupUng 
is stronger at infinity than the single particle potential. If this is not satisfied, our proof does not 
work. There could be some physical reason behind this. If a stationary state exists, this means that 
even if the chain is in a state of very high energy, the mean time to reach a region with low energy 
is finite (see e.g. [Has80]). But if m < n, the relative strength of the coupling versus the one-body 
potential goes to zero at high energy. The consequence is that there is almost no energy transmitted 
between particles. Since the only points where dissipation occurs are the ends of the chain, we see 
that the higher the energy of the chain is, the slower this energy will be dissipated. Probably this is 
not sufficient to destroy the existence of a stationary state, but it could explain why the proof does 
not work in this situation. It is even possible that this phenomenon destroys the compactness of 
the resolvent of K. 



7 The invariant measure 

This section is devoted to the proof of Theorem 2.7. Throughout this section, we denote by T* 
the semi-group generated by the system of stochastic differential equations (2.7). We also assume 
that A1-A3 are satisfied, so Proposition 2.5 and 2.6 hold, as well as Theorem 2.4. The proof 
of Theorem 2.7 is divided into three separate propositions, showing respectively the following 
properties of the invariant measure /x: 

(i) Existence and smoothness. 

(ii) Decay properties. 

(iii) Uniqueness and strict positivity. 

Proposition 7.1. If Assumptions A1-A3 are satisfied, the Markov process given by (2. 7) possesses 
an invariant measure /i. It has a density h, which is a C°° function on R^^+^. 

Proof. By Theorem 2.4, we know that K has compact resolvent. This implies also the compact- 
ness of the resolvent of and thus of Lq. Since G grows algebraically at infinity, we see that 
the constant function 1 belongs to Hq. Moreover, we notice that LqI = 0, thus the operator Lq 
has an eigenvalue 0, which is isolated because of the compactness of its resolvent. This in turn 
implies that Lg also has an isolated eigenvalue 0. We denote the corresponding eigenvector by g 
and normalize it so that {g, 1)^^ = 1. Since Lq is hypoelliptic, g must be C°°. 
Assume first that g>Q. We then define 



where Zq is the normalization constant appearing in the definition of Ti^. Set ji{dx) = h{x) dx; we 
want to check that fi is the invariant measure we are looking for. Notice that //(dx) is a probability 
measure because 



h{p,q,r) = Zq ^g{p,q,r)e 



-2/3„G 



(7.1) 




1 . 
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Let A be a Borel set of R^^+^. Then the characteristic function xa of ^ belongs to Hq. We have 

((T*)>)(yl) = J {T\^){x)fi{dx) = Z,'J e-^PoOi.^)g{x){T\^){x)dx 
= Z^' j e-^^oGi^){{T^yg){x)xA{x) dx = /.(A) , 



thus /X is an invariant measure for the Markov process defined by (2.7). 

The argument showing that it was indeed justified to assume g positive can be taken over from 
[EPR99a, Prop. 3.6]. □ 

We next turn to the decay properties of the invariant measure h. We first introduce a convenient 
family of Hilbert spaces. 

Definition 7.2. Choose 7 e R. Wfe define the Hilbert space W'^'^'> as 

= \3{X, (a;) dx) = V{G^) . 
We will denote by (■, ■)(-y) cind \\ ■ ||(^) the corresponding scalar product and norm. We also define 

7>0 

which is the set of all functions that decay at infinity faster than any polynomial. 

We already know that his a C°° function, so we want to show that it is possible to write 

h{p, q, r) = h{p, q, r)e-^oG(p,9,0 ^ ^ |y(oo) _ 

The function h being an eigenfunction of the operator K, the decay properties of the invariant 
measure are a consequence of the following result. 

Proposition 7.3. The eigenfunctions of K and K* belong to C°°{X) fl W^°°\ 

We will show Proposition 7.3 only for the eigenfunctions of K. It is a simple exercise left 
to the reader to retrace the proof for the eigenfunctions of K*. We already know that K and K* 
are hypoelliptic, so their eigenfunctions belong to C°°{X). It remains to be proven that they also 
belong to 

To prove the proposition, we will show the implication 

/eiyw and x/eiyw feW^^+'K (7.2) 

which immediately implies that the eigenvectors of K belong to W'^°°\ For this purpose, we 
introduce the family of operators defined by 

: V{K^) ^ 

where T>{K^) is given by 

V{K^) = {/ G I Kf e W^W} . 

The expression Kf has to be understood in the sense of distributions. 
We have the following preliminary result. 

Lemma 7.4. C^{X) is a core for K^. 

Proof. The proof uses the tools developed in Appendix B and is postponed to Appendix C. □ 
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The invariant measure 



The key lemma for the proof of Proposition 7.3 is the following. 

Lemma 7.5. There are an e > and constants > such that for every 7 > and every 
u e V{K^), the relation 

holds. 

Proof. Since we know that C^{X) is a core for K^, it suffices to show (7.3) for u G C^{X). Let 
L be the first-order differential operator associated to a divergence-free vector field. Then we have 

for/,^eCo°°, 

(L/,^)(^) = -U^LG^^g) = -{f,G'^Lg) - 2^{f,G'^G-\LG)g) 
= -(/,^^)(,)-27(/,G-i(LG)^)(,). 

We write this symbolically as 

L; ^ -L^ - 2-fG-\LG) . 
We can use the latter equality to show that there are constants c^^ and c^^ such that 



2 ^"^"^ 



Using the explicit form of K, this in turn yields the useful relation 



(7) -^L\\^r^"'\\{-y) ^R\\'-'rj^<^\\(^) 

(7) + (^RlW'R - ^R'iNm\(^) -r \'«, Jx"/(7) 



+ 411 (^L - ^l(1o)u\\L) + 411 (r^ - \RqN)u\\L) + {uJku), 



where ff^ is some bounded function. 

We now have the tools to prove the validity of (7.3). We use Proposition 2.5 to write 



(7.4) 



2^_^^^ = WG'G-fuW^ < CiWKG-^uW^ + ||G"^ii|n 
< C{\\Gmuf + \\[K, G-^M^ + \\G^u\\^) . 

An explicit computation yields 

[K, G'y]u = G^ {h d,^ + \ + /o)« , 
for some smooth bounded functions and /q. We are thus able to write 

||«||?,+,) < C{\\Ku\\l^ + \\u\\l^ + cm^uWl^ + cm.^^Wh) ■ (7.5) 
Using (7.4), we can write 

<C{\\Ku\\l^ + \\u\\l^). 
This, together with (7.5), completes the proof of the assertion. □ 
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Proof of Proposition 7.3. Lemma 7.5 immediately shows that T>{K^) c W^'^'^^^ for every 7 > 0. 
This proves the assertion (7.2). 

Let / be an eigenfunction of K. We know that / G L^(A') and, because it is an eigenvector of 
K, we have Kf e Ll Thus, by (7.2), / e W^'\ Of course Kf e PF^") as well, so / e W^^'\ 
This can be continued ad infinitum, and so we have / e W'^°°\ which is the desired result. □ 

Finally, we want to show the strict positivity and the uniqueness of the invariant measure. 
The proof of this result will only be sketched, as it simply retraces the proof of Theorem 3.6 in 
[EPR99b]. 

Proposition 7.6. The density h of the invariant measure /i is a strictly positive function. Moreover, 
the invariant measure is unique. 

Sketch of proof. The idea is to show that the control system associated with the stochastic differ- 
ential equation (2.7) is strongly completely controllable. This means that, given an initial condition 
Xq, a time r and an endpoint x^, it is possible to find a realization of the Wiener process w such that 
^(r; Xq, w) = x^. The main assumption needed to show that is that the gradient of the two-body 
potential is a diffeomorphism. This is ensured by assumption A3. 

The consequence is that, for every time r, every initial condition Xq and every open set U, the 
transition probability P(t, Xq, U) is strictly positive. Because ji is invariant, we have 

li{U) = j P{t, X, U) n{dx) > . 

This implies the strict positivity of h. Uniqueness follows from an elementary ergodicity argument. 
□ 



A Proof of Lemma 4.4 

Throughout this appendix, we will make use of the same notations as in Section 4, i.e. Ti. = 1^ (^'*)> 
V = C^(R") and D is the set of differential operators with smooth coefficients. 

Moreover, A denotes some finite subset of D and is identified with closed operators on H. The 
operator is defined as 

A2 = 1 + J^AM. (A.l) 

AeA 

We will moreover assume that HI and H3 concerning A and ^ holds, i.e. A,BeA and / G 5" 
imply 

[A,B]ey\A), A*ey^A), [AJ]ed. (A.2) 
In order to prepare the proof of Lemma 4.4, we need a few auxiliary results. 

Lemma A.l. Let A, 'S, T> and A be as above and assume HI and H3 hold. Then, if A e ^^(.4), 
the operator AA~^ is bounded. 
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The proof of this lemma will be a consequence of 

Lemma A.2. Let A, ^, V and A be as above and assume HI and H3 hold. Then, ifAi, A2 G A, 
the operators A^A'^ and AiA2A~^ are bounded. 

Proof. Let us show first that A^A~^ is bounded. Since "D is a core for A, it suffices to show that 
there is a constant C such that 

\\Ajr<c\\Afr v/ep. 

This is an immediate consequence of 

iiA/ii^ = ii/ip+$]p/ir 

In order to show that AiA2A~'^ is bounded, we will show that there are constants r and C such that 

K^,/|P<C||AV + (r-l)/ir (A.3) 
We can write the following equality 

\\(A'-l)f + rfr^rW + 2rJ2\\Afr+ E {f,A*AB*Bf) 

AgA A,BeA 

= rW + 2r \\Afr + iW^Bff + (/, [A* A, B*]Bf)) . 

AeA A,BeA 

We can write the operator intervening in the last term as 

[A* A, B*]B = A*[A, B*]B + [B, A]*AB . 
Because of HI, this implies that there are positive constants Cj^^q such that 

\\{A' - i)f +rfr >rw +2TY,wr + E 11^^/11' 

AgA B,CeA 

- E C^Bc\\Af\\\\BCf\\. 

A,B,CeA 



If we use now 

s 

we see that we can choose r big enough to have 



2xy < x^s^ + ^ x,y >0 , s > 



AeA B,CeA 

This immediately implies (A.3). □ 
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This lemma can now be used to prove Lemma A. 1 . 

Proof of Lemma A.l. We want to show that A e ^^(.4) implies AA~^ bounded. We already 
treated the cases i = 1 and i = 2. For the other cases, we proceed by induction. Let us fix j > 2 
and assume the assertion has been proved for i < j. Then the operators of the form 

^1^2 • ... • AjA-^ A^eA, (A.4) 

are bounded. We distinguish two cases. 

j = 2n. We write the operator of (A.4) as 

A,A,A-' ■ A^A,A,A-^ • . . . • A'-' A,^_,A,^A-'- . 
We show that operators of the form 

^2m-2^^^-2m A, B E A , m<n, 

are bounded. We write 

The first term is bounded by Lemma A.2. The second term is bounded by noticing that 
[\'^rn-2 ^ AB] e yf^'^(A) and using the induction hypothesis. 

j = 2n + 1. We write the operator of (A.4) as 

A^A^A-^ ■ A^Ag^A-^ • . . . • A2"A2n+i^"^''"^ ■ 

The first terms are bounded exactly the same way as before. Concerning the last term, we 
have 



^2ji+1^ ^ 



[A2", A. 



2n+l 



which is bounded by Lemma A.2 and the induction hypothesis, noticing that the commu- 
tator belongs to ^^^"(.A). 

This completes the proof of the lemma. □ 

We need another result from [EPR99a]. 

Lemma A.3. Let {A{z)} C ^(7i) be a family of uniformly bounded operators, A > 1 a self- 
adjoint operator and let F{X, z) be a real, positive bounded function. Then 



f 

Jo 



A{z)F{A,z)fdz 



< sup P(|/) II 11/11 

y>0 



supF(A, z) dz , 

A>1 



(A.5) 



If furthermore A — A{z) is independent of z, one has the bound 



AF{A,z)fdz 



< 



sup 

A>1 Jo 



F{X,z)dz, yfen. 



(A.6) 
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Lemma A.4. Let A, and A be as above and assume HI and H3 hold. If X & 3^5(«^)' then the 
operators 

A'^XA^ with l3 + -f <-j 

are bounded. 

IfY&£is such that [Y, A^] e y^{A), then the operators 

A^[A'*,y]AT' with a + P + -f <2- j 

are bounded. 

IfX,Y G 2, are such that 

[X, A'] e y^{A) , [Y, A'] e yl{A) and [[A', X],Y] e yi-^'-'iA) , 

then the operators 

A^[[A°,A:],y]AT with a + (3 + < A - j - k 

are bounded. 

Proof. Let us prove the first assertion. The case 7 = is handled by noticing that 

A(^X = A^+^{X*A-^y , 

and that both operators of the latter product are bounded by Lemma A. 1 . The case P — Ois handled 
in the same way by considering the adjoint. 

The proof for the other cases follows exactly [EPR99a]. We will demonstrate the techniques 
involved by proving the third assertion, assuming the first two assertions hold. The second assertion 
can be proved in a similar way without using the third one. 

We will first assume that a e (—2, 0). In this case, we can write (see e.g. [KatSO, § V.3.11]) 

A- = r z-l\z + A2)-i dz, = _!!^W2) _ 

We notice moreover that it is possible to write 

%z + A^Y\X\Y\ = + A^)-^ [[A^ X\Y\ {z + A^)-^ 

+ + A2)-i[A2, X\{z + A2)-i[A2, Y\{z + A2)-i (A.8) 
+ + A2)-i[A2, Y\{z + A2)-i[A2, X\{z + A2)-i . 

If we substitute the expression (A.7) in A'^ [[A", X\ , F] A''^ and use (A.8), we get three terms, which 
we call T^, and T3, and which will be estimated separately. 

Term T^. This term is given by 
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We define B = [[A^, X^Y] e y^+^-'^{A) and write 



AT 



/■oo 



oo \p A 7 

[A2 , S] . dz 



z 



+ A2' 



(^ + A2)2 



The term T^^ is estimated by writing, for any f G H, 



\TiJ\\ 



Jo 



[z + A^y 



2-j-k-l3\ 







sup 




A>1 . 


/O 




roo 


sup 




A>1 . 


10 



dz 



(Z + A2)2 



Since the assumption yields B e y^+'^ the norm is bounded. The integral is also bounded 

because, by assumption, we have a + ^ + P<4: — j — k. 

To bound T^2' we observe that [A^, B] e y^+''-^{A). Using (A.6), we find the bound 



W2. 



A^ 



< 



sup 

y>o 



Z + A2 

A^ 



y + A^ 



[A2, B]A^- 
[A2,5]A3- 



-j-k-13 



-fdz 



-j-k^f3 



(^ + A2)2 

oo ^^l+P+j+k-S 

z"/^ sup 







A>1 (z + x^Y 



dz 



This expression is bounded when a + P + ^<4:— j — k and a e (—2,0). This can be seen by 
making as before the substitution z A^s. 

Before we go on, we introduce the notation A^ = {z + A^)"^. 

Term T2. This term is given by 



T2 = C, 



z""/^ —B — dz 



z + A^ z + A^ z + A^ 
where we defined 

A=[A^,X] and B = [A^,Y]. 
Since [A^., B] — A^[B, A^jA^, the term appearing under the integral can be written as 

A^A^AA^SA^AT = A^A^ASA^AT + A^A^AAJS, A^JA^A^ . 
According to this, the term is split into two terms T21 and 722- We have 



|7^2i/ll < 11/11 sup \\Af'A^ABA-f'-^-'^\ 

y>0 



roo 

/ s"/2 sup 

Jo A>1 



{s + iy 



ds . 



The integral is bounded by hypothesis. The norm is also bounded, because AB e y^~^''{A). For 
the second term, we have 



1^22/11 < 11/11 sup \\A^AyAA^[A',B]A-^-^-'\ 
y>o 



^o+/3+7+j+fc-4 
5"/^ sup : — — TTT^ ds . 



\>1 



{s + iy 
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This is bounded in the same fashion, noticing that 

sup \\K^KyAKy[K^,B]K-^-i-''\\ < sup ||A^A^AA-^-^|| sup 

Term T3. It can be bounded in the same way as by symmetry. 

We now have to check the assertion for the other values of a. If o; = or o; = 2, it holds 
trivially. For a > 0, we proceed by induction, using the equality 

[[A"+2, X],Y]^ h? [[A", X], y] + A" [[A2, X], r] + [A°, y][A2, X] + [A2, y][A", X] . (A.9) 

For a = —2, the assertion is proved using equality (A. 8) with z = Q. For a < —2, we also proceed 
by induction, using (A.9) with 2 replaced by —2. This completes the proof of Lemma 4.4. □ 



-A^^+/^-2[A2,S]A-' 



B Proof of Proposition 1.2 

Proposition B.l. T*, as defined in (1.10), extends uniquely to a quasi-bounded strongly continu- 
ous semi-group on dx). Its generator L acts like C on functions in C^{X). 

Proof. See the proof of Lemma A. 1 in [EPR99a]. □ 

We now turn to the question of the domain of the generator L. Recall that £ is the set of all 
formal expressions of the form 

^ai{x)D\ k>0, aeC°°(R"). 

\l\<k 

To any element L e SI having the above form, we associate its formal adjoint L* e £ in an obvious 
way. In the sequel, the notation (/, g) will be used to denote the scalar product in L^ if f,g G L^ 
and the evaluation f{g) if / is a distribution and g e C(f (R"). We hope this slight ambiguity will 
not be too misleading. 

We associate to every L e £ the operator : V{Tj^) by 

{TLf){x) = Lf{x) and P(rj = {/ e L^ | L/ e L^} , 

where Lf has to be understood in the sense of distributions, i.e. 

{Lf){g) = f{L*g) for all geC^{R-). 

We also define the operator Sl : T>{Sl) L2(R") by 

The operators and are usually called the minimal operator and the maximal operator con- 
structed from the formal operator L. The following result is classical, so we do not give its proof 
here 
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Proposition B.2. For every L E Z, we have Tl — S^* and — Tj^». In particular, this shows 
that T]^ is closed. □ 



We prove now the quasi m-dissipativity of -S'^. We define 

M 



i=l 



By definition, if Sji is strictly m-dissipative, is quasi m-dissipative. It is well-known that an 
equivalent characterization of strict m-dissipativity is that 

(a) is strictly dissipative and 

(b) Range(5^) = H. 

Proposition B.3. Assume AO holds. Then is strictly m-dissipative. 

Remark. It is clear that the statement holds if we consider the minimal operator in L^(/C, dx), 
where /C is some compact domain of X. The idea is to approximate A" by a sequence of increasing 
compact domains and to control the rest terms. 

This proposition fills a gap in [EPR99a], since the statement "Re(/, L*f) = — |||cr^V/||^ -|- 
(/, div h f) < in the proof of Lemma A.l is not justified for every / e V{L*). 



Proof. Property (a) is immediate. By the closed-range theorem, property (b) is equivalent to the 
statement 

(b') / e and C*f = imply / = 0. 

Assume on the contrary that there exists a non-vanishing 
function / e for which £*/ = holds in the sense of dis- 
tributions. Since C* is hypoelliptic, / must be a C°° function. 
Let us choose some function 93 e C^(R_,_) such that (p{x) = 1 
if X e [0, 1]. We also define 



X 




R 

^{G{x)/n) . 
By assumption, C* f — 0, so we have 

= 2Re(^„/, tf) = {^J, tf) + (£*/, ^nf) ■ 

Since e Co° the other functions are C°°, we can make all the formal manipulations we 

want. In particular, we have 



</, {^^t + C^^)f) = . 



(B.l) 



Recall that £ is given by 

MM MM 



1=1 

M 



i=l 



1 . 



1=1 



i=l 



(B.2) 



i=l 



56 



Proof of Proposition 1.2 



Straightforward computation yields 

M M 



1=1 i=l 
M M 



(B.3) 

M 

1=1 ' 

M 



i=l 

Using AO, we next verify that | < C* for all a; G A" and for all n > 1. We define 

Ci = supcp"{x) and C2 = supxcp'{x) . 

x>0 x>0 

An elementary computation shows that AO implies that there exist constants Cg, . . . ,05 > for 
which 

\d^G{x) \ < C3 , \d^G{x)\'^ < c^G{x) , and (r^ - \jFi{p,q)f < c^G{p,q,r) . 
We thus have 



i=l * 
/ C1C3 + C2C4 7iC2C5 \ ^ ^ 



=1 

as asserted. Moreover, the first part of AO implies that there exist constants C, a > such that 

supp$,„ d {xeX \ ||a;||" > n/C} . (B.4) 
Substituting (B.3) back into (B.l), we get 

M 



= -2f^C||y^5../|r-|lv^/||'+ / 



Since f &l3{X), one has 

|2 _ II ^112 



(B.5) 



n— >oo ' ■ 

Moreover, the uniform boundedness of together with property (B.4) imply that 

lim / $„(a;)|/(a;)|2dx = 0. 

This supplies the required contradiction to (B.5), thus establishing the strict m-dissipativity of S^. 

□ 
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We complete now the 

Proof of Proposition 1.2. It only remains to be proved that L — and that L* — S^* . 

It is clear that the generator L of T* satisfies C L. Since -S'^ is quasi m-dissipative, i.e. has 
no proper quasi dissipative extension, and since the generator of a quasi-bounded semi-group is 
always quasi m-dissipative, we must have L = Sjr. 

Concerning the adjoint, we have by Proposition B.2, L* = T^*. It is possible to retrace the 
above argument for C* to show that S'^. is quasi m-dissipative. Since L* is also quasi m-dissipative 
and S'^* C L*, we must have L* — Sj^*. □ 



C Proof of Lemma 7.4 

Using the technique developed in Appendix B, we can now turn to the proof of Lemma 7.4. Recall 
that K is given by (2. 10) and that 

Moreover, is the maximal operator constructed from K when considering it as a differential 
operator m M^W. We have 

Proposition C.l. Cq'{X) is a core for K^. 

Proof. We introduce the unitary operator U : W^'^'^ — > L^(A') defined by 

{Uf){x)^G^{x)f{x). 



We also define = K^, \ C^{X). The operators K^^ and are unitarily equivalent to the 
operators K and respectively by the following relations. 



b 



u 






u 















K.. 



U~ 



u 





u ^ u 















7 U 



By construction, is maximal. Thus, by Proposition B.2, its adjoint K* is minimal. It is imme- 
diate that the formal expressions for K* and are given by 

K; = G-m*G^ and = G^KG-^ . 

It is now a simple exercise to retrace the proof of Proposition B.3 to see that K* and are both 
m-accretive. The remark of Section 1.2 concerning the adjoints of m-accretive operators implies 
that is also m-accretive. Since C K^, we must have — and thus — K^. This 
proves the assertion. □ 
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